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Linearized fluctuations of quantized matter fields and the spacetime geometry around 
de Sitter space are considered in the case that the matter fields are conformally invariant. 
Taking the unperturbed state of the matter to be the de Sitter invariant Bunch-Davies state, 
the linear variation of the stress tensor about its self-consistent mean value serves as a source 
for fluctuations in the geometry through the semi-classical Einstein equations. This linear 
response framework is used to investigate both the importance of quantum backreaction 
and the validity of the semi-classical approximation in cosmology. The full variation of the 
stress tensor 5(T a b ) contains two kinds of terms: (1) those that depend explicitly upon the 
linearized metric variation 5g c d through the ([T a b ,T cd ]) causal response function; and (2) 
state dependent variations, independent of 8g cc i- For perturbations of the first kind, the 
criterion for the validity of the semi-classical approximation in de Sitter space is satisfied for 
fluctuations on all scales well below the Planck scale. The perturbations of the second kind 
contain additional massless scalar degrees of freedom associated with changes of state of the 
fields on the cosmological horizon scale. These scalar degrees of freedom arise necessarily 



from the local auxiliary field form of the effective action associated with the trace anomaly, 



are potentially large on the horizon scale, and therefore can lead to substantial non-linear 
quantum backreaction effects in cosmology. 



PACS numbers: 04.62,+v, 95.36.+X, 98.80.Qc 



I. INTRODUCTION 



Although matter is undeniably quantum, in General Relativity its gravitational effects are treated 
classically. In order to include the quantum effects of matter, a natural next step from the purely 
classical theory is a semi-classical treatment, in which the classical matter stress tensor T a b in Einstein's 
equations is replaced by its expectation or average value {T a b ), but the spacetime geometry is still 
treated classically. This is clearly an approximation to a more complete treatment, whose validity 
depends upon a number of assumptions, chief among them that the quantum fluctuations of the matter 
energy-momentum-stress tensor T a b about its average value are "small," at least on macroscopic distance 
scales. 

This heuristic condition can be given a precise meaning in the semi-classical or mean field limit of 
an effective field theory approach to gravity. In the limit, ^ — ^ but the number of quantum fields 
N — > oo with hN fixed, it becomes permissible to replace the conserved T a b of a quantum theory by its 
expectation value (T a ,), as the source term for the semi-classical Einstein equations, 



The corrections to the renormalized expectation value {T a b ) R are suppressed by 1/N, and negligible in 
the large N limit [1^ . In this way and in this limit the treatment of the metric geometry of spacetime 
as classical, sharply peaked about its mean value, can be obtained from an underlying quantum theory 
of matter. 

Replacing the classical source of Einstein's equations by a quantum expectation value has a number 
of important implications. Since T"i{x) is an operator which is at least bi-linear in local quantum fields 
at the same spacetime point, any direct evaluation of its expectation value is divergent. The expectation 
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value (T a b ) must be renormalized, i.e. its divergent parts identified and removed by means of suitable 
local counterterms in the effective action of low energy gravity. 

It is essential that the renormalization procedure maintain the covariant conservation of (T a b ) R , for 



it is to act as a source of the gravitational field through the semi-classical Einstein equations (1.1); 



otherwise (1.1) would be inconsistent. Covariant methods for identifying and removing the ultraviolet 
divergences {T a b ) by point-splitting, heat kernel expansions, or dimensional regularization which satisfy 
the requirements of general covariance have been developed, and lead to consistent well-defined finite 
results for {T a b ) R [21|3|. 



Eqs. (1.1) also provide consistency conditions for the finite terms in the renormalized expectation 
value {T a b ) R . To evaluate this quantity in the usual Minkowski vacuum state of flat spacetime, one is 
required to define the renormalized expectation value of T a b to vanish identically, when A = 0, or (1.1 ) 
will not be satisfied for R a b = 0. Equivalently, by transposing the cosmological term A to the right 



side of (1.1), the existence of a flat spacetime solution to the semi-classical Einstein equations requires 
that the difference, (T a b ) R — 5 a b vanish for R a b = 0, so that any cosmological term is exactly 



canceled by a finite renormalization of {T a b ) R in flat space. Thus, neither the vacuum expectation value 
(T a b ) R nor the cosmological term A 5 a b alone have direct physical significance, but only the combination, 
(T a b) R ~ sJg ^6) which acts as the source of curvature. 



The condition that flat spacetime be a solution to the semi-classical Eqs. (1.1) with all matter in 
its Lorentz invariant Minkowski vacuum state is a renormalization condition on the low energy effective 
theory of gravity. This condition serves to define the effective theory by fixing the subtractions of 
infinite and finite vacuum terms, and makes the semi-classical framework a predictive one in spaces 
other than flat space and/or in states other than the vacuum state. At present this renormalization 
condition cannot be justified by any fundamental theory, but solely by appeal to the experimental fact 
that the very high frequency fluctuations responsible for the divergences in (T a b ) do not produce any 
appreciable mean curvature of space. The semi-classical theory described by ( |1.1[ ) must be regarded 
then as a low energy effective theory [3] , on a par with pion effective theories of nuclear forces or the 
Fermi theory of the weak interactions, whose UV divergences are subsumed into effective low energy 
parameters, and whose predictions are reliable only at sufficiently low energies or long wavelengths, 
pending a more complete, but unspecified theory at high energies and short wavelengths. 

Classically, with A > and in the absence of any other sources, there is a maximally symmetric non- 
flat solution to Einstein's equations, namely de Sitter spacetime with R a b = A 5 a b . In the classical theory 
A is a fixed constant with absolutely no dynamics, and de Sitter space is a stable solution to Einstein's 



3 



equations with a positive cosmological constant. On the other hand, in the semi-classical framework 



of (1.1), as just observed, the cosmological constant term cannot be divorced from the definition of 
{T a b ) R itself, as A can be reabsorbed into a finite redefinition of {T a b ) R . Since {T a b ) R depends upon the 
state in which it is evaluated, the energy of the vacuum is no longer a fixed non-dynamical constant of 
the Lagrangian, as in the classical theory, but instead depends upon the infrared choice of "vacuum" 
state in curved space, in which the expectation value (TV) R is defined. Which infrared state and which 
spacetime is preferred then becomes a dynamical question, which can be investigated by probing the 
response of the system to perturbations of the state of the fields and the spacetime geometry together 



on length scales much greater than the Planck scale. In the semi-classical theory (1.1) the stability of 



de Sitter space to these infrared perturbations depends upon the behavior of the quantum fluctuations 



of the matter sources on cosmological length scales of the order of H 1 = -y/3/A, and the dynamical 



question can be studied in a systematic way with the same effective action that gives rise to (1.1) by 
the method of linear response. 

The linear response of a system to a small perturbation is familiar from other branches of physics, 
and has been formulated for semi-classical gravity in Refs. [3 [6]. The fluctuations of the matter fields 
lead to fluctuations in the expectation value (TV) R , which is probed by considering the response of 
{T a b ) R to a small perturbation of the geometry g a b — > g a f, + 5g a b- The quantum action principle tells us 
that this variation is given formally by 

S(T%(x)) = \J d A x'y^gJ^)U a b cd(mt) (x,x') 5gM) 

= -y/ dV7^M^t,t')(m|[n(x),r cd (x')]|zn) 6 9cd (x') , (1-2) 

up to contact terms proportional to delta functions and derivatives thereof with support at x = x' 
[6]. The local contact terms are determined by the renormalization of the retarded polarization 
IL a b (x, x') of the matter field(s), evaluated in the same state as the self-consistent background 



solution of ( 1.1 ). 



The linear variation (1.2) is essentially equivalent to first order perturbation theory in the metric 
perturbation 5g a b = h a b, and assumes that the entire change in the quantum state and the expectation 
value, (T a b ) comes about directly from the change in the geometry. These variations, linear in the 
metric variation h a t>, we term variations of the first kind. This first class of linear perturbations are the 
ones upon which attention has been focused in the earlier literature on this subject [TJ El El [TUl ITT1 [T2] . 
However it is also possible to consider variations of the state that lead to variations in (TV) R which 
are not directly driven by local variations in the metric. For example one might consider variations of 



the temperature of a system coupled to a heat reservoir, or other variations of the boundary conditions 
imposed on the system. These variations may depend upon dynamical degrees of freedom contributing 
to the stress tensor other than the local metric geometry, and hence are not proportional to Sg c d- The 
corresponding state dependent variations in (T a b ) R which do not depend upon the local Sg c d, we term 
variations of the second kind. We shall see that the trace anomaly parameterizes certain specific state 
dependent variations of the stress tensor which are of this second kind, and that they can have significant 
effects in de Sitter space. 



The dynamical linear response equation which probes the self-consistent solution of (1.1) to small 
fluctuations of both kinds is the linear integro-differential equation for the self-consistent metric pertur- 
bation hab, 

S J fl « _ | n + A5 a\ = 87tGn5 {T a )a ? (L3) 



obtained when the local linearized variation of the left side of (1.1 ) is set equal to the properly renor- 
malized right side. Once the ultraviolet divergences of this linear response equation are isolated and 



removed, (1.3) can be used to study the infrared fluctuations of matter and geometry of de Sitter space 
in a self-consistent and reliable way. This is our main purpose in this paper. 



The most direct way of deriving (1.3) and determining the correct renormalization counterterms 
necessary to properly define the formal expressions for (T a .) and 6(T°l) is to integrate out the matter 
fields, obtaining their one-loop effective action in a general gravitational background. In the large 
iV limit, the contributions of higher order gravitational metric fluctuations themselves to 5{T a b ) are 
suppressed by 1/N and may be neglected at leading order. In this limit the metric may be treated as a 
(semi-) classical quantity. In particular, there are no stochastic or "noise" terms in the linear response 
equations, ( |1.3| ) at leading order in the 1/N expansion. These can arise only at higher orders and are 
outside the domain of the semi-classical theory. 

Renormalizing the one-loop matter effective action and taking its first variation yields the semi- 



classical Einstein equations (1.1). In addition to fixing the background geometry, the solution of (1.1) 
involves also fixing the quantum state (or density matrix) \in){in\ of the fields, with respect to which 
all expectation values are to be computed. Then a second variation of the effective action, or first 



variation of ( 1.1 ) is performed around this background metric to yield the linear response equation ( 1.3 ) 



for h a b- The kernel of the polarization integral in ( 1.2 ) is computed using the background geometry and 
the Green's functions for the quantum field in the specific quantum state (density matrix) fixed in the 
previous step. The form of the counterterms used to define the renormalized expectation value {T a b ) R 
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are the same used to renormalize its linear variation 5{T a b ) R |6j. 



In [6] we proposed a criterion for the validity of the semi-classical Einstein equations (1.1 ), namely, 



that the solutions of the linearized equations (1.3) should lead to linearized gauge invariant amplitudes 
which remain bounded for all times. The criterion is satisfied in flat spacetime. The present work 
extends the analysis and tests the criterion in de Sitter spacetime for conformal matter, where the 



variation of the stress tensor in (1.3) also can be evaluated in closed form. 

The interest of applying this linear response method to de Sitter spacetime arises from several differ- 
ent considerations. First de Sitter spacetime serves as a fundamental testing ground for semi-classical 
effects of quantum vacuum fluctuations in a curved space with a cosmological constant, in which analytic 
results can be obtained. Secondly, current cosmological models of structure formation in the universe 
assume that this structure grew by gravitational instability from small quantum inhomogeneities in 
an otherwise featureless, primordial de Sitter-like inflationary phase. Finally, observations of Type la 
supernovae suggest that the universe may be in a de Sitter-like accelerating phase in the present epoch, 
with some 70-75% of the energy density of the universe today being of a hitherto undetected dark variety, 
with negative pressure, p ~ — p [13] . Since this is exactly the equation of state of a cosmological term 
or that of the quantum "vacuum," the study of quantum fluctuations and their stability in de Sitter 
spacetime is of possible direct relevance to the cosmology of the early universe, structure formation, 
and the cosmic acceleration of the present universe. 



In this paper we study the solutions of this dynamical linear response equation (1.3) under small 
fluctuations in de Sitter spacetime for the case of conformal matter fields, i.e. those which are classically 
conformally invariant, up to the quantum conformal or trace anomaly. The study of metric perturba- 
tions in cosmology due to quantum conformal matter in cosmology was initiated in Ref. [7], whose 
authors made use of an in-out effective action formalism, rather than the retarded boundary conditions 
appropriate for causal linear response. The in-in formalism for the effective action in cosmology was 
developed in Ref. [11] , and extended and elaborated in Ref. [12 . The general form of the perturbed 
stress tensor in an arbitrary FRW cosmological spacetime obtained in |12| . agrees with the earlier re- 



sults of Refs. [H [10]. We study the solution of the resulting causal linear response Eqs. (1.3) with 
this general stress tensor source for inhomogeneous and anisotropic perturbations away from de Sitter 
space in Sec. VII A. Previous discussions of linearized perturbations of de Sitter space have focused 
instead on spatially homogeneous, isotropic perturbations, earlier in the functional Schrodinger initial 
state formalism |14j . and more recently for a minimally coupled scalar field in the in-in effective action 
formalism |15 . 
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For the more general case of spatially inhomogeneous and anisotropic perturbations, we find the 
analogous short distance (i.e. Planck scale) quantum fluctuations expected from the corresponding 
analysis in flat spacetime, which should be discarded as outside the range of validity of the semi- 
classical theory. Going further, when general inhomogeneous and anisotropic state dependent variations 
of the stress tensor are considered (i.e. variations of the second kind), we find additional modes arising 
from coherent macroscopic changes of state on every scale, characterized by a spatial wave vector k, 
including the scale of the de Sitter horizon itself. These cosmological horizon modes are associated with 
the effective action of the trace anomaly, not variations of the local metric geometry, and as such are 
different than anything in the purely classical theory, or encountered in previous studies of restricted 
classes of perturbations. The existence of this new set of scalar cosmological horizon modes arising from 
the quantum trace anomaly is a principal result of this paper. 

To keep this result in perspective, and when considering quantum effects in cosmology generally, it is 
important to recognize that the intrinsically quantum phenomenon of phase coherence and consequently 
non-classical effects may be present at any scale. Indeed macroscopic quantum states are encountered 
at low enough temperatures in virtually all branches of physics on a very wide variety of scales. In 
addition to low temperature and condensed matter laboratory systems, the chiral symmetry breaking 
expectation value of quark bilinears (qq) in QCD and the Higgs field ($) in electroweak theory are 
examples of macroscopic vacuum condensates, which extend over arbitrarily large distances. It is these 
condensates which provide the connection between a more microscopic approach to both the strong and 
electroweak interactions and the low energy effective theories which preceded them. Such coherence 
effects due to the quantum wave-like properties of matter, and its propensity to form phase correlated 
states over macroscopic distance scales cannot be treated in a purely classical description of gravity. 
However, the semi-classical treatment of the effective stress-energy tensor source {T a b ) R for Einstein's 
equations already allows for such effects, in a mean field description |16j . 

A significant part of the effective action for conformally invariant fields is determined by the conformal 
or trace anomaly. Although intrinsically non-local in terms of the spacetime geometry, the effective 
action determined by the anomaly can be cast into a covariant local form by the introduction of one or 
more scalar auxiliary degree(s) of freedom [161 El (IS HHJ [201 [21] . Since there are two distinct cocycles in 
the non-trivial cohomology of the Weyl group in four dimensions [22 , the most general representation 
of the anomaly action is in terms of two auxiliary scalar degrees of freedom, each satisfying fourth 



order linear differential equations of motion (4.8). The resulting stress tensor in terms of these scalar 
auxiliary fields has been used to study the non-local macroscopic quantum coherence effects contained 
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in the semi-classical effective theory (1.1) |16| [23] . The scalar auxiliary fields are related to the freedom 
to change the macroscopic quantum "vacuum " state, and supply additional infrared relevant terms in 
the low energy effective theory of gravity [22J. In perturbations about de Sitter spacetime we shall show 
they give rise to new cosmological horizon modes, not contained in classical perturbation theory, that 
can lead to large backreaction effects at the very largest Hubble scale of a de Sitter universe. Thus scalar 
degrees of freedom in cosmology arise naturally from the effective action of the conformal anomaly of 
massless quantum fields in the Standard Model, without the ad hoc introduction of an inflaton field. 

The paper is organized as follows. In order to fix ideas and notations we review in the next section 
the linear response fluctuation analysis in flat spacetime with a vanishing cosmological constant. In 
Section III, we consider linear response for conformal matter fields in conformally flat spacetimes, 



where the stress tensor variation of the first kind (1.2) can be expressed in closed form. We show 
that for long wavelength solutions of the equations within the range of validity of the semi-classical 
theory, it is sufficient to study only the time-time component of the linear response equations, in a 
general Friedman-Robertson- Walker (FRW) spacetime. In Section IV we consider the stress tensor and 
linear response equations that follow from the trace anomaly part of the effective action only, in the 
auxiliary field description, comparing the result to the previous general approach in conformally flat 
spacetimes. In Section V we derive the linear response equations in both approaches for fluctuations 
in de Sitter spacetime, and show how the auxiliary field degrees of freedom in the anomaly action 
describe cosmological horizon scale modes. In Section VI we discuss coordinate transformations and 
recast our linear response equations in a completely gauge invariant form, showing in particular that 
the new cosmological horizon modes are gauge invariant. We give the quadratic gauge invariant action 
corresponding to the linear response equations in these variables. In Section VII we give a qualitative and 
quantitative discussion of the solutions of the linear response equations in de Sitter space, showing that 
the auxiliary field anomaly action captures all the infrared quantum effects correctly. In Section VIII the 
cosmological horizon modes, their interpretation as fluctuations in the Hawking-de Sitter temperature, 
and their possible implications for the the stability of de Sitter space are discussed. Section IX contains 
a full summary and discussion of our results. There are two appendices, the first containing an in-depth 
discussion of the properties of the non-local kernel K which arises in the linear response equations in 
conformally flat spaces, and the second, some of the details of the second variation of the anomaly 
action. Except where indicated the units used are h = c = 1 and the metric and curvature conventions 
are the same as those of Misner, Thorne, and Wheeler 
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II. LINEAR RESPONSE IN FLAT SPACETIME 



In this section the derivation and solutions to the linear response equations in a flat space background 
are reviewed. In flat spacetime Lorentz invariance permits the full decomposition of perturbations 
into scalar and tensor components, with respect to the background Minkowski four-metric r] ab = diag 
( — h ++). In Ref. [6] explicit results for these components were obtained for a quantum scalar field 
with general mass m and curvature coupling £. Making use of those results as a specific example, we 
generalize the discussion here to any set of underlying quantum fields. We emphasize that the four 
dimensional tensor decomposition in Minkowski spacetime is distinct from the decomposition of the 
perturbations into scalar, vector, and tensor components with respect to the spatial three-metric (/in- 
appropriate for general FRW spacetimes, and employed in the succeeding sections. 

The starting point of a linear response analysis of quantum fluctuations is a self-consistent solution 



to the semi-classical Einstein equations, (1.1). In flat spacetime this requires that (T a b ) R = in the 
usual Minkowski vacuum state with A = 0. Since it is a dimension four operator, the renormalization of 
(T a b ) requires fourth order counterterms in the effective action. These gives rise to additional geometric 



terms in equations (1.1), fourth order in derivatives, which are usually displayed on the left side of 



the semi-classical equations (1.1). We shall instead adopt the convention of including such geometric 



terms in the definition of the renormalized expectation value, (T a b ) R itself: cf. equations (2.1) and 



(2.2) below. Viewing the semi-classical theory as an effective field theory, these local higher derivative 
terms multiplied by finite coefficients are in any case negligibly small at macroscopic length scales much 
greater than the extreme microscopic Planck scale, Lpi = \J KG at / c 3 ~ 1.616 x 10 -33 cm., or at energy 
scales much less than the Planck energy M Pl c 2 = ^/hc 5 /G N ~ 1.223 x 10 19 GeV. 



In flat spacetime, the non-local integro-differential equation ( 1.2 ) in coordinate space is easily handled 



by Fourier transforming to momentum space. With the self-consistent solution of (1.1) for A = just 



the usual Minkowski vacuum state, the linear response equations around flat spacetime take the form [6] 

5G ab = 87TG N 5(T ab ) R 

= 8nG N | (a R + 5 ^H ab + + 5 W# a6 j . (2.1) 

Here the two local, conserved geometric tensors, ^H ab and ^H ab , also denoted as — A ab and — B ab 



9 



respectively are 

^H ab = ^ J d 4 x C abcd C abcd = 4V c V d C (a c h f + 2C a \ d R cd = -A ab , (2.2a) 
1 S f 

^H ab =—— b J d 4 x^gR 2 = 2g ab UR-2V a V b R + 2RR ab - g -^R 2 = -B ab . (2.2b) 



These tensors, which are derived from the two dimension-four local invariants, C abcc iC abcd and R? 
respectively, are required for renormalization of the expectation value {T ab ). The parameters a R , and 
(3 R are the finite, renormalized coefficients of these local terms. The causal response functions, 

F ™ ( K^K*r% , (2.3) 

Jo s s + — te sgnuj 



in (2.1) have been renormalized by subtracting the dispersion integral three times at s = 0. 1 They 
are the non-local part of 5{T ab ) R arising from the matter fluctuations. The renormalization effected 
by subtracting the Taylor series expansion in K 2 of the unrenormalized response function up to order 
(K 2 ) 2 fixes the value of the cosmological constant term, Newtonian constant (order K 2 ) and fourth 
order terms A ab and B ab at order (K 2 ) 2 to be their prescribed renormalized values at K 2 = 0. The 
renormalization prescription and local counterterms of the retarded polarization tensor U a b cd ^ ret \x,x') 
at x = x 1 are defined in momentum space by this procedure at K 2 = 0. The modification of the 
renormalization procedure necessary for massless quantum fields is discussed below. 



The superscripts T and S in (2.3 ) denote quantities associated with tensor (spin-2) and scalar (spin-0) 



perturbations which are defined with respect to the four dimensional background Minkowski metric by 



Eqs. (2.14) below. As discussed in Ref. [6] these causal response functions are completely determined 
by their imaginary parts in terms of the spectral functions p( T ' S ^ of the matter fields. The spectral 
functions p^ T ' S ^ are finite and may be computed directly from the one-loop cut diagram of Fig. [l]in flat 



spacetime. 



It is essential that the linear response equation (2.1) is causal, i.e. that the stress tensor response at 



(t, x) is sensitive to variations of the geometry only within the past light cone of (t, x). The — ie sgnw 
prescription (where K 2 = — lo 2 + k 2 , k = \k\ and u = K°) enforces retarded boundary in-in conditions, 
appropriate for this causal linear response. 



Note also that both terms on the right side of (2.1) are geometrical, involving the linearized metric 



fluctuation h ab itself, i.e. they are homogeneous variations of the first kind only. It is clear that 

1 In this section only we employ the notation K a for a four-vector in flat spacetime, and K 2 = r/ a bK a K b — —lo 2 + k 2 to 
distinguish it from k = \k\, the magnitude of its three spatial components used throughout the paper. 
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FIG. 1: The one-loop vacuum polarization H al f d (x, x'). Its imaginary part, obtained by cutting the internal lines 
as illustrated and placing them on-shell gives the spectral functions p( T,s >. 



non-geometrical terms could always be added to 5(T a b ) R , by allowing states other than the Minkowski 
vacuum to contribute to the variation, even in the absence of any metric variation. These state dependent 
contributions will be considered in the next section. 



The expressions (2.1)-(2.3) are finite and general for the linear response of a quantum field theory 



in the usual Minkowski vacuum perturbed around flat spacetime. To linear order all indices are raised 
and lowered with the flat space metric 77^. The dependence on the particular quantum matter theory 
enters only through the spectral functions p^ T ' S \ which must be positive definite on general grounds of 
unitarity. Explicit calculations of the spectral functions and causal response functions of a free scalar 
field theory with mass m and curvature coupling £ are given in Ref. [B]. The results are: 

„n 2 

(2.4a) 



P {S) (s) 



6(s - Am 2 



24vr 2 
3 — Am 2 



1 



4m 2 



m H 



1 



Am 2 



m 

A 



(2.4b) 



60vr 2 V s 

For conformal scalars, m = and £ = g, the scalar spectral function p^ vanishes for all s > 0. However, 
the expression in (2.3) is finite in the limit m — > limit. The two response functions in this limit 



arc 



F (S) 
F (T) 



1440vr 2 

1 

960vr 2 



In 



Ab 
K 2 



le sgn uj 



2b In 



K 2 



te sgn uj 



(2.5a) 
(2.5b) 



where b 



(i^p" 120 



is the coefficient of the C a bcdC abcd term of the trace anomaly for the conformal scalar 



field in an arbitrary curved background. The logarithmic divergence in (2.5b) in the limit of zero mass 



is the result of the renormalization of a R at K 2 = 0. If the renormalization is performed instead at an 
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arbitrary but non-zero mass scale fx 2 , and we define 



a R (/j, 2 ) = a R + bin 



m- 



then the combinations appearing in the linear response equation (2.1) are 

f( t ) 



12 



a R (fj> ) + b In 



- le sgnu; 
~J 2 



(2.6) 

(2.7a) 
(2.7b) 



These are finite in the conformal limit m — > 0, provided a R (fj, 2 ) remains finite in that limit. Note also 
that since 

K 2 



ft 



d 
dfj, 2 



or(m ) + bln 



re sgn to 



0, 



(2. 



the physical linear response does not depend on the arbitrary renormalization scale /i 2 , i.e. this combi- 
nation is renormalization group invariant. For the scalar response function of a conformal matter 
field, there is no logarithmic term and no fj? dependence at one-loop order, and the last term in (2.7b) 



is simply an additional finite renormalization of the coefficient j3 R of the R 2 term in the effective action. 



Eqs. (2.5 ) and (2.7) are valid in fact for any conformal matter field, provided one uses the appropriate 



anomaly b coefficient for that field. This may be seen from the general form of the trace anomaly for 
conformal matter fields in an arbitrary curved spacetime, viz., 



{T a a ) R = bF + b' \E--UR \ + b"UR 

b[F+ 2 -UR\+b>E+(b" _ 26 ± 26 ") DR . 



where 



771 — *r> *jyabcd p jyabcd a tj jyab , r>2 

= tXabcd tt — H a bcd-tl — 4/t a {,/t + H 



and 



t~i z-iabcd 

= ^abcd^ 



p jjabcd n p r>ab . 



R 2 



(2.9) 

(2.10a) 
(2.10b) 



with Rabcd the Riemann curvature tensor, *R a bcd = \ e abefR e ^ cl i its dual, and C a b c d the Weyl conformal 
tensor. The coefficients b, b' , and b" are computed at one-loop order, proportional to h, and determined 
by the number of massless conformal fields and their spin via 

n 



b 

y 



120(4tt) 2 

h 

~360(4vr) 



(N s + 6N F + 12A^ y ) , 
(N s + UN F + 62 AV) , 



(2.11a) 
(2.11b) 
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with Ng the number of spin fields, Np the number of spin \ Dirac fields, and Ny the number of 



spin 1 fields [3_. The h in (2.11) emphasizes that these coefficients appear at one-loop order, and that 
the anomaly terms are present in the leading order semi-classical limit, h — * 0, N — > do, with hN finite 
(where N is any of N$, Np, Ny). With this understood, we shall henceforth set H = 1. 

The b coefficient is also the one that determines the logarithmic scale dependence of the C abc dC abcd 
term in the effective action. The coefficient b" is scheme dependent, corresponding as it does to the trace 
of the tensor ^H ab = —Bab, which is derived from a local action, whose coefficient can be shifted by a 
finite shift of (5 R . For classically conformal invariant theories there is no logarithmic scale dependence 
in the coefficient of the R 2 term in the effective action at one-loop order and one finds b" = 2(6 + b')/3 
in covariant regularization schemes, such as dimensional regularization, so that the last term in ( |2.9| ) 
vanishes 125 . 



From these considerations we deduce that the renormalized stress tensor variation in (2.1 ) in Fourier 
space for conformal matter in flat spacetime may be written in the form, 



S ( T ab) R , co . 



nf. 



ot R (n 2 ) + b In 



K 2 



te sgn oj 



5A 



ab 



(2.12) 



which is obtained by substituting Eq. (2.7) into Eq. (2.1). Since both the invariants E and F are 
quadratic in the curvature tensor, their linear variations away from flat spacetime vanish, and only the 



!&d<5.R term survives in the variation of the trace, 



S(T a a ) 



R, conf. 



2 b 



6/3 fl + y ) n(SR). 



(2.13) 



All of the dependence of the linear response on the conformal matter theory in flat spacetime is contained 



therefore in the single parameter b, determined by the trace anomaly (2.9) and given by (2.11a) for free 
conformal matter fields of any spin. The b' anomaly coefficient of the Euler- Gauss-Bonnet density E 
does not enter the linear variation of the stress tensor around flat spacetime. 



The tensor and scalar components in (2.1) or (2.12) are conveniently separated with the aid of the 
set of orthogonal projectors, 

1 



P. 



P 



AT) cd 
ab 

(V)cd 
ab 

(S) cd 
ab 

Jab 



] a c e b d + e a % c 



i 

2 



')cd 



1ab^ 



■yd 



Vab -d a —d b 



flab 



KgK b 

K 2 



(2.14a) 
(2.14b) 
(2.14c) 
(2.14d) 
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in momentum space, so that the linearized tensor and scalar metric perturbations around flat spacetime 
are 



h (T,S) _ p (T,S)cd h 
h ab = F nh ri cd , 



ab 



(2.15) 



respectively. Both perturbations are transverse in the four dimensional sense, 



b h (T,S) 



d b h 



ab 



= K°h 



b h (T,S) 



ab ' 



(2.16) 



while the tensor perturbation is also tracefree in the four dimensional sense, 



abAT) 



rTh 



ab 



0. 



(2.17) 



The linear variations around flat spacetime of the three local geometric tensors appearing in (2.1) may 
be expressed then as 

1 



5G ab 
5A ab 
SB ab 



*d£ + SG ab 



2 U ' L ab + U a ab 



2 ab 



K 2 h (S) 

A n ab ' 



m 2 h[ s b ^ 6(**) 2 *2\ 



(2.18a) 
(2.18b) 
(2.18c) 



in terms of this decomposition into tensor and scalar metric components in momentum space. The 



remaining projector onto linearized vector perturbations (2.14b) is not needed because h 



i iV) 
l ab 



P. 



(V)cd 



ab 



cd 



is gauge dependent and non-transverse, and hence cannot appear in the linear variations of the conserved 



tensors in (2.1 ). 



Making use of the projection operators (2.14), we find that the linear response equation around flat 



spacetime (2.1) separates finally into two independent, orthogonal components, viz., 

K 2 



8G 



5G 



(T) 
ab 

(S) 
ab 



h 



(T) 
ab 



8ttG n < a R {n z ) + b In 



1 



:0 ab 8R 



-8nG N (6/3 fl + |) (K 2 ) 2 h[ S b \ 



ab ' 



(2.19a) 
(2.19b) 



The solutions of the linear response equations in flat spacetime were discussed in ref. [6] for scalar 
fields with arbitrary mass and curvature coupling £. In Fourier space it is straightforward to show that 



there are no low energy solutions to (2.19) other than the usual propagating transverse, traceless, spin-2 
gravitational wave modes at K 2 = of the pure Einstein theory with no sources. Since these transverse 
degrees of freedom are massless, there are only two helicity states which propagate, the other three 



spin-2 states and the one apparent scalar mode in (2.19b) being eliminated by the four diffeomorphism 
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constraints coming from the ti and tt components of the equations in a canonical analysis splitting 
spacetime into space + time. 

Additional non-trivial solutions to both the semi-classical tensor and scalar equations appear at K 2 ~ 
Mpj (for a R and /3 R of order unity). However, these solutions are outside the range of applicability of the 
semi-classical effective theory, which requires K 2 <C M^,,, and hence should be discarded. Physically, 
this must be the case since there is no independent length scale against which the wavelength of the 
metric perturbations of a conformally invariant quantum field can be compared in flat spacetime other 
than the Planck length. The situation will be quite different for cosmological spacetimes where the 



horizon scale plays an important role. From (2.13) or (2.19b) it is clear that the Planck scale solutions 
in the scalar or trace sector involve 5R 7^ 0. These Planck scale solutions will be present in the linear 
response equations around arbitrary spacetimes as well. Hence we shall find it convenient to set 5R = 0, 
in order to exclude these Planck scale solutions in the effective field theory description, and focus on 
the remaining solutions (if any) that are determined by the low energy curvature or horizon scale in 
cosmology, which is assumed many orders of magnitude greater than the Planck scale. 



III. LINEAR RESPONSE OF CONFORMAL MATTER IN COSMOLOGY 

In this section some general properties of conformally invariant fields in conformally flat spacetimes 
are briefly reviewed along with some previous computations of the perturbed stress-energy tensor for 
these fields in these spacetimes. Then the specific form of the linear response equations in a general 
spatially flat Friedman-Robertson- Walker (FRW) cosmological spacetime is derived and the infrared 
physical content of the equations is shown to reside purely in the tt component. 

A conformally flat spacetime is one for which the metric is 

g a b{x) = VL 2 {x)ri ab (3.1) 

with rjai, the Minkowski metric. The function Q(x) is called the conformal factor. Conformally invariant 
fields are ones whose classical action is invariant under conformal transformations when the field is 
multiplied by some power of the conformal factor. If the spacetime under consideration is conformally 
flat and the matter fields are conformally invariant, then the state of the fields that is most often 
considered is called the conformal "vacuum" [3i. It is obtained by mapping the Minkowski vacuum 
and all of its Green's functions in flat spacetime to the conformally flat spacetime using the same 
transformation of the fields which keeps the action invariant under conformal transformations. In this 
special state the expectation value of the stress tensor of the conformally invariant fields is completely 
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determined by its trace anomaly (2.9) and can be written in terms of local curvature invariants and 



their derivatives in the geometry specified by (3.1) in the form [3) 126] . 

b 



(Tab) 



R ,conf 



-2b' ®H ab --8 fl6 



(3.2) 



where 



ab 



2 l 1 2 

-R^R c b + - R R a b + - R c dR c 9ab — -^R g ab — 2 C ac M R 



cd 



(3.3) 



and B a b is defined by (2.2b). Of course, just as there are many states in flat space that are of interest 



besides the Minkowski vacuum, the conformal state is simply one choice of allowed state among many, 
even if we restrict ourselves to the subclass of states which are spatially homogeneous and isotropic, 



consistent with the symmetries of (3.1). Applying the term, "vacuum" to the conformal state is also 



somewhat misleading, since {T a i > ) R conf / 0. 



Some time ago several authors considered the linear variation of the expectation value (3.2) for 



conformal matter fields by conformally transforming the linear variation (2.12) in flat space j9j[T0]. It 
was assumed in that work that the fields are in the conformal vacuum state and the variation considered 
was of the first kind, namely it was assumed that the quantum state of the fields follows the metric 
variation in a prescribed way, with no other state dependent variations considered. The results of Refs. 
[21 [TU] are used in this section to write the specific form of the linear response equations in a spatially 
flat FRW spacetime, with the addition of the state dependent terms which are omitted in [91 IIP] , but 
certainly allowed on general grounds [3]. 

The line elements for spatially flat Friedman-Robertson- Walker (FRW) spacetimes can be written in 
the alternate forms, 



ds 2 



-dt 2 + a 2 {t) dx 2 = n 2 (r]) (-drj 2 + dx 2 ) 



(3.4) 



where 7] is the conformal time and t is the comoving time. These cosmological spacetimes are conformally 
flat with the conformal factor £l(r]) = a(t), the FRW scale factor. 

In the flat space version of the linear response equations it was useful to perform a full Fourier 
transform. For a FRW spacetime it is useful to Fourier transform only in space, but not in time. The 



inverse Fourier transform can be used to undo the time part of the logarithmic kernel in (2.12 ) with the 
result, 



K{j] -rj';k;n) 



dio 
2^ 



-co 2 + k 2 



it sgn u> 



(3.5) 
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This is a distribution, which may be defined by its action upon integration against a suitably smooth, 
restricted class of test functions f(rf) for which J dr]'K(rj — Tf)f(rf) exists. Because of the — iesgnu; 
prescription, K has support only for times rf earlier than the time 77, consistent with causal linear 



response. At rj = rf the distribution (3.5) is singular and requires regularization. The properties of this 



distribution and two possible regularizations are studied in detail in Appendix A. 

In terms of this distribution the linearized perturbation of {T ab ) R for a conformal field in the conformal 



vacuum around a conformally flat spacetime with metric (3.4) is (9j ITO] 



8(T ab (ir,k)}j 



b_ 
02 



drj'K(rj - rf\ k; /j 2 ) SA ab (rj'; k) +L ab + S(T ab (r]; k)) R . 



(3.6) 



Here L nb is the local contribution, 



L a b 



-a R {^) 5A ab -[p R + l] 8B ab - 2b' 5 <%, 



(lnfl) 5C a 



c d 



(3.7) 



and the final term in ( 3.6 ) is the variation of the stress tensor arising from the variation of the state away 



from the conformal "vacuum," which was omitted in j9j [TO], but which is present in the more general 



context of state dependent variations [3 . Except for the lower limit of the integral in (3.6) introduced 
to allow for the possibility of starting the linear response system at some finite initial conformal time 



770, and the addition of the variation of the final state dependent term in (3.6), these are the results for 



the variation of the stress tensor of conformal matter around a conformally flat space as given by the 
authors of Refs. [SHU]. 



In (3.6) both 5A ab and SC ' c > are evaluated on the linearized perturbation, 



h ab = 2 h ab . 



(3.8) 



from flat spacetime, with h ab the full linearized metric perturbation about the conformally flat spacetime 
g ab = £l 2 i~i ab of (3.4). The linearized variation around flat space of ^ c >H ab = — A ab is given in momentum 
space by Eqs. (2.14)-(2.18) of the previous section, while 5C c b is given by the completely traceless 
part of 

1 



SR. 



c d 



d b d c h a d + d a d d h b c - d a d b h ca - d c d d h ab 



cd 



3C ad 1 



(3.9) 



in the usual way, all derivatives and indices referring to flat Minkowski coordinates. 

As in the corresponding expression for flat spacetime (2.12), the a R (fi 2 )5A ab term in (3.7) comes 



from varying the local C a bcdP Weyl invariant term in the effective action, while the f3 R 5B ab term is 
a finite renormalization of the local R 2 term in the effective action, which we have chosen to include in 
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the variation of {T a b) R . The b 5B a b and b' 5 ^ 3 « a 6 terms are the variations of the background geometric 
terms in (3.2), while the d c d^\hinSC „ c 6 d ] term has been obtained by the authors of Refs. [9] [10] in 



different ways, and is necessary for consistency. This term will be shown in the next section to be 
a necessary consequence of the covariant effective action functional determined by the trace anomaly 



(2.9). 



Since SA a b, 5A a b and the last term in (3.7) involving the Weyl tensor variation are traceless, the 
trace of the linear response equations for conformal matter perturbed away from a FRW background 
for which C a b c d = simplifies considerably. In fact, for variations of the first kind, the trace is deter- 



mined completely by the trace anomaly (2.9), which is given in terms of local curvature terms and is 
independent of the quantum state of the field. Hence the trace of the variation is given by the local 
geometric equation, 

- SR= 8irG N 8(T a a ) R = 16nG N j-26' (^R a c 5R c a - ^RSRj + \3(3 R + ^ <5(D-R)| . (3.10) 
By assumption, in the semi-classical approximation we require the Ricci curvature and its variations 



to be much smaller than the Planck scale. Thus all the terms on the right side of (3.10) are small 
compared to the left side. Hence we may restrict ourselves to only those solutions of the full linear 
response equations for which 

5R = 0. (3.11) 

The only solutions eliminated by this condition are Planck scale solutions, such as the ones we have 
already encountered in flat spacetime, with k 2 ~ Mj, ; , which lie outside the range of validity of the 
semi-classical approximation under consideration. Henceforth we will restrict ourselves therefore to the 



tracefree or non-conformal solutions of the linear response Eqs. (1.3) around FRW spacetimes. Earlier 
authors have studied the conformal variations of the metric and stress tensor for conformally invariant 
fields in de Sitter space in more detail, and found no interesting gauge invariant modes other than those 



expected on the Planck scale [T3]. The condition (3.11) eliminates these trace or conformal Planck 
scale solutions from the start, and simplifies the analysis of the remaining tracefree, non-conformal 
perturbations. This condition also eliminates the class of quantum inflationary solutions considered in 
Refs. |27l [28]. 

With 5R = 0, we are restricted to the sector of tracefree metric perturbations (in the four dimensional 
sense). Further, in expanding about a FRW background one can decompose the metric perturbations 
into scalar, vector and tensor perturbations with respect to the conformally flat three dimensional 
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metric g%j = a 2 7]ij, and focus on the scalar sector (in the three dimensional sense). The vector and 
tensor perturbations for conformal matter around a conformally flat spacetime are not expected to 
show any non-trivial degrees of freedom in linear response other than those in the standard classical 
theory. 

The metric perturbations which are scalar with respect to the background three-metric can be 
parameterized in terms of four functions, (A,B,C,£), of the form \29\ [30] . 



ht 
h tj 
hij 



-2A 



adjB — * iakjB 



2a z 



\ 3 



(3.12a) 
(3.12b) 
(3.12c) 



in momentum space. As is well known, and elaborated further in Sec. VI, only two linear combinations 



of these four functions are gauge invariant. Since we have also required (3.11), this means that there 



remains only one dynamical gauge invariant metric function to be determined by linear response in this 
scalar sector. The information about this remaining metric degree of freedom is contained completely 
in the tt component of the full linear response equations. 



To show this, let the background semi-classical Einstein Eqs. (1.1 ) be written in the form, 

R 



n = <n>« 



i 



R a 



b~ ^ 5a b + A 



0. 



(3.13) 



Their covariant conservation, 



V7 q-a _ f) n-a _i_ pa q-c _ -pc q-a _ n. 
V a J- h — O a l b -\- 1 ac l b 1 b l c — U , 



(3.14) 



and first variation imply 



d t {5T\) + ft(*Tj) + P jt {5T\) - r 3t {5T\) = , 

d t {5T\) + d^STi) + T^ t {5T\) - P lt (5T)) - r* y = . 



(3.15a) 
(3.15b) 



for the time and space components respectively. Here the background equations (3.13) have been used 



along with the spatial homogeneity and isotropy of the background. This leaves the Christoffel symbols 



shown in (3.15) as the only remaining non-zero ones. In coordinates (3.4) the non-vanishing Christoffel 



symbols are 



1 jt - 1 tj 



ij ~ a ■na 



(3.16a) 
(3.16b) 
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where the overdot denotes differentiation with respect to the comoving proper time variable t. Note that 



because of the self-consistent background equations (3.13), one may equally well consider the variations, 
STi, 5T a i) or 5T ab , freely raising and lowering indices with the background metric. 



Using relations (3.16) in (3.15a), and imposing the time-time component of the linear response 
equations, 



sr\ = o 



(3.17) 



and the total trace of the linear response equations, 



6T a a = 5T\ + 5T\ = 0, 



(3.18) 



equivalent to (3.10), gives 



di(6T\) = 0. 



(3.19) 



For scalar perturbations 5T\ = g l3 dfP for some scalar function V, there being no way for a spatial vector 



to appear in the scalar sector other than by differentiation with respect to x % . Then (3.19) implies that 
X7 2 V = or V = up to a possible constant. Since a constant in V is eliminated in 5T\ = g^djV, 



(3.19) implies 



sr\ = o, 



(3.20) 



or in other words, the time-space components of the linear response equations around a general FRW 



space are automatically satisfied for scalar metric perturbations of the form (3.12), provided the self- 



consistent background equation (3.13), the time-time component (3.17), and the trace equation (3.18), 
and are all satisfied. 



To complete the proof one then notes from (3.17) and (3.18), that the spatial trace 5T\ vanishes, 



leaving only a possible traceless part, in 5T l - which for scalar perturbations can be expressed in the 
form, 



ST) 



1. u„ 1 



(3.21) 



for some scalar function Q. Substituting this form into (3.15b), and using (3.17), and (3.20) gives 



d%Q = 0, which by an argument identical to that for the time-space component implies 



5T) = 0. 



(3.22) 
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Thus, (3.20) and (3.22) show that the sufficient conditions for all of the semi-classical linear response 



equations to be satisfied around a general FRW background for scalar perturbations of the metric (3.12 ) 



arc: 



(i) the background is self-consistent, so that Eq. (3.13) is satisfied; 



(ii) the time-time component (3.17) is imposed 



(iii) the trace equation (3.18) is imposed 



Since the trace equation can be imposed by (3.10), for non-Planck scale fluctuations, the tt component 



contains the only remaining gauge invariant dynamical information about scalar perturbations parame- 



terized by (3.12) around a self-consistent solution of the semi-classical Einstein equations (3.13). Being 



able to focus solely on the time-time component will simplify our analysis considerably. 

Since the variation of (T a b) of conformal matter in a conformally flat spacetime is determined by the 



form of the trace anomaly (2.9), another route to the linear response equations in this case is afforded 



by the effective action and stress tensor corresponding to the anomaly, which has been derived in several 
previous works. In the next section this alternative derivation from the effective action and stress tensor 
of the anomaly is given and compared with the above derivation based on Refs. [9j [TO] . 



IV. ANOMALY ACTION AND STRESS TENSOR 



For massless matter or radiation fields, there is no intrinsic length scale associated with their quantum 
fluctuations, and the conformal variation of the one-loop effective action in a general curved background 
is determined completely by the trace anomaly. Although the anomaly determines uniquely only the 
part of the effective action which responds to conformal variations, and leaves the remaining confor- 
mally invariant part undetermined, it has been argued elsewhere that the anomaly action contains the 
only effective field theory corrections to Einstein's theory which can be relevant in the infrared, i.e. 
at macroscopic distance scales much greater than the Planck scale j22[ [3~T] . The availability of the 
alternative formulation of linear response for conformal matter around conformally flat backgrounds 
reviewed in the previous section allows this argument to be tested in FRW spacetimes by comparing 
directly the energy-momentum tensor and linear response equations derived from the anomaly action 
with the previous exact formulation. 

The general form of the effective action which records the effects of the trace anomaly is non-local 
in the metric. However in several earlier works a local form was obtained by the introduction of two 



21 



scalar auxiliary fields, tp and -0 |16| [23] . This auxiliary field effective action is of the form, 

Sanom = b' S d ^ om [g; ip] + b S'^mb; V> V>] , (4.1) 



with 



S^L\g-M = l J d 4 x^~g L(n V ) 2 + 2(R ab -^yV a ^(V b <p)+(E-^DR\ <p\ (4.2a) 
S£ljg;<p,i>}= J d^x^-g {-(□^(□V0+2(V-|^) (V a ^)(V 6 ^) 



+ \f V +\(e-IuR\^\ . (4.2b) 



2 r 2 V 3 

By varying this coordinate invariant effective action functional with respect to the spacetime metric, 
we obtain the stress-energy tensor, 

T ~ = b'E ab + bF ab (4.3) 

where the two separately conserved tensors are give by Eqs. (3.41) and (3.42) of |16j . namely, 

E ab = -2(V (a <p)(V b) n<p) + 2V c [(V c ^)(V a V^)] - ^V a V fe (V(^) 2 
+ \ Rab (V<^) 2 - 4 R c ia {V b) <p){V c p) + \R (V a <p)(Vb<p) 
+ ^g ab {-3(D</?) 2 + □ {Vtp) 2 ) + 2 (3R cd - Rg cd ) (V c¥ >)(V^)} 

c d ' 



- - V a V fe □ ip - 4 C a c b a V c V d tp - 4 R c {a V b) V c <p +-R ab Up+-R V a V b <p 
2 

~3 



{V {a R)V b) ip+^g ab {2n 2 ip + 6R cd V c V d ip-4Rnp + (V c R)V c p} , (4.4) 



and 



F ab = -2 (V (a ^)(V 6) n*l>) - 2 (V (o ^)(V 6) U<p) + 2 V c [(V c ^)(V a V fe V) + (V c V)(V a V^)] 
-\ V a V 6 [(V c ^)(V c V0] + ^ Rab (V c ^)(V c ?/0 - 4 i? c (a [(V 6) p)(V c V) + (V 6 )^)(V C ^)] 



+^i2(v (a¥? )(v b) v) + ^5ab{ -3(n^)(nv)+ □ [(v c ^)(v c ^)] 

+2 (3i? cd - i? 5 crf ) (V c ^)(V d 0)} - 4 V c V d (C (0 C 6) V) - 2 C a c ^i?, 



■ | v a v 6 □ ^ - 4 c a c 6 d v c v d v - 4 i? c (o V 6) V c + ^i?o6n^ + -i? v a v 6 v 
-|(V (a fl)V 6) ^+^«6 {2n 2 + 6 J R cd V c V^-4i?nV + (V c J R)(V c V)} , (4.5) 



where (V^) 2 = (V ' a ip){V a p). By varying (14. 1|) with respect to the local auxiliary fields one obtains 
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their equations of motion, 



A 4 ^ 



E 

~2 
F 

2~ 



UR 

~3~ 
1 



CabcdC 



abed 



(4.6a) 
(4.6b) 



which are linear in ip and ip. Here 



A 4 = + 2R ab V a V b - ^RU + ^{V a R)V a ■ 



(4.7) 



Since the terms quadratic in the auxiliary fields in (4.2) are conformally invariant, they lead to only 



traceless terms in the corresponding stress tensors (4.4) and (4.5). Therefore the terms in the trace are 



linear in ip and ip, and in fact are proportional to (4.6), yielding the local geometrical traces 



E a a = 2A A <p = E- -OR, 
F a a = 2A 4 V = F = C abcd C 



abed 



(4.8a) 
(4.8b) 



Thus the classical trace of the stress tensor following from the local effective action (4.2) reproduces the 



quantum trace anomaly (2.9). 



Moreover, it was shown in Refs. |21[ [22] [32] by considering its conformal variation under g ab — > g ab 
e~ 2a g ab that 



-g $g 



ab 



T wz [g;a]=2m ab+ 1 -WH ab , 



where T^ z [g; a] is the stress tensor of the Wess-Zumino effective action, 



v wz WM = I 



d xyf— ~g 



2aA 4 a+ [ E- -BR ) a 



(4.9) 



(4.10) 



obtained by varying the metric g ab , keeping a fixed. Considering equations (4.1)-(4.4), dropping the 



overbar on (4.9) and using the notation of (2.2) one finds that 

E ab 



\ip=2 In n 

in the full metric g ab , provided 



L ab 



9\° 



2 ^H ab + -B ab , 



9ab = e v r] ab = £l 2 r]ab 



(4.11) 



(4.12) 



is conformally flat, with the auxiliary field ip = 21nS7. It is easily checked from the conformal variation, 



E- -UR = 4A 4 lnft = 2A 4 <^ 
3 



(4.13) 
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in the full metric, that ip = 21nQ satisfies (4.8a). These relations may be proven directly as well, by 



the use of formulae for the conformal variations of the Ricci tensor, such as 



Rab 

R = 



-VaV b (f - i ( V aV ?) (VfeV?) + ^ [- □ V9 + ( Vp) 



(4.14a) 
(4.14b) 



in the full metric g a b 

Thus, provided that the auxiliary field is fixed in terms of the geometric conformal factor via tp 



21nf2, the tensor E ab reduces to the purely local combination of geometric tensors given by (4.11). 
Taking the second auxiliary field tp = 0, and the Weyl tensor C a c b d = 0, gives 



rpanom 
I ab 



b'E ab = -2b' ®H ab + b - B ab 



(4.15) 



for the full contribution of the anomaly effective action (4.1) to the expectation value of the energy- 
momentum tensor of conformal matter in an arbitrary conformally flat spacetime. We note that for 
a conformally flat spacetime of the form (3.4) and ip(t) = 2 In a(t) + Co, the tensor E a b derived from 
the auxiliary field effective action gives rise to precisely the geometric tensor ^H a b discussed in Ref. 
which is not derivable from a local geometric action and for that reason was called "accidentally 



conserved" in [3]. The actual origin of this tensor is the form of the anomaly action (4.2). 



Eq. (4.15) differs from the known exact result (3.2) only by the addition of a b"\3R term to the 
anomaly and corresponding local R 2 term to the effective action with b" = 2(6 + b')/3, the same choice 



of b" expected from the flat space result: cf. Eqs. (2.9)-(2.13). Since such local terms with arbitrary 



coefficients should be added to the effective action in any case, we conclude that the action (4.1) plus 
these additional local terms gives the correct general form for the stress tensor of conformal matter in 
its conformal "vacuum" state in an arbitrary conformally flat spacetime. 



To study linear response next we wish to vary the anomaly stress tensor (4.3) away from the confor- 



mally flat background. Since the relations (4.11 ) and (4.15 ) hold for arbitrary conformal transformations, 
it follows that they hold also for linearized conformal transformations about a conformally flat FRW 



spacetime. Hence it is clear that in the variation of (4.15) we must expect to obtain local geometri- 



cal terms which are just the variations of the local geometrical terms in (4.15). By adding the local 
F = C a hcdC abcd and R 2 terms to the effective action, it is clear that one should expect the variation of 
the local geometrical tensors A ab and B ab which arise from them. Thus the origin of all the local terms 



in (3.6) and (3.7) is clear, except perhaps for the last term in (3.7), whose origin we now discuss. 
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Inspection of the cp dependent terms in (4.5) shows that 



SF ab 



-8V c V d 



(lnJ2) 5C {a ^ -4(lnn)dC {a c b fR cd 



(4.16) 



around any conformally fiat spacetime with C c b d = 0. Since the Lagrangian density ^/—gFip from 
which this term is derived is conformally invariant under g a b = £l 2 r] a b - > Vab, ^(^ab) anc ^ ^he cov ariant 



derivatives in (4.16) may be replaced by their values, <5C/ C fe ^ and d c dd on the perturbation (3.12) in flat 



spacetime by a simple overall scaling with n 2 . That is, one may check explicitly that in general 



2V c V d 



V C (ab) + <P C (ab) R cd 



(a b) 



1 

n 2 



2V c V d 



-~ c d 



(a b) Red} 



(4.17) 



for any two metrics conformally related by g a b = Q 2 g a b- Varying this relation keeping 93 and n fixed, 
with 5g a b = Q 2 5g a b, and taking g a b = r\ a b to be flat gives upon making use of (4.16), 



5F, 



ab 



92=21nQ,i/>=0 



n 2 



d c dd 



(lnn)C (a ' 



(a b) 



(4.18) 



Hence the last local term in (3.7) of the variation of the stress tensor expectation value of conformal 



matter in a conformally flat spacetime may be regarded as a direct result of the anomaly effective action 



(4.1)-(4.2) in the auxiliary field formalism. Since the A ab tensor is just proportional to (4.17) with ip 



replaced by a constant, we also deduce that 



n 4 



A a , 

A b 1 



(4.19) 



upon raising one index, for the two metrics related by (4.12). 



Thus the anomaly action and stress tensor reproduces all the local terms denoted by L a b in (3.7) in 



the exact result (3.6). However it does not reproduce the non-local term in the variation (3.6), involving 



the logarithmic distribution K. This term comes from the part of the one-loop effective action which is 
invariant under conformal transformations g a b — > Q 2 g a b: and does not contribute to the anomaly or the 
anomaly action. It may be argued from the addition of the local a R (fi 2 )C a b c dC abcd term, renormalized 
at an arbitrary mass scale [x that a non-local logarithmic term involving K with the correct coefficient 
must be present in the effective action in order for the scale /i 2 to drop out of the low energy dynamics 



according to (2.8). Whether the non-local term is essential or not in perturbations around de Sitter 



space will be studied in detail in Sec. VII. 

Through the anomaly effective action and its auxiliary field description, it is possible to express the 
variation of the stress tensor in a completely local form, and to gain the advantage of examining a wider 
class of state variations parameterized by the auxiliary fields as well. By specializing to de Sitter space 
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in the next section it is shown that the independent variation of the metric and auxiliary fields of the 
anomaly stress tensor (4.3) exhibits specific additional state dependent horizon scale degrees of freedom 
which are not present in the classical Einstein theory. 



V. LINEAR RESPONSE OF CONFORMAL MATTER IN DE SITTER SPACETIME 

In this section we specialize the discussion to de Sitter spacetime, and derive the linear response equa- 
tions in the anomaly action auxiliary field formulation, comparing it in detail to the general expression 



for the linearized stress tensor variation in (3.6). 



In FRW coordinates with flat spatial sections (3.4) the scale factor of de Sitter spacetime is 



a(t) =e Ht = -^ r = n( V ) , T) € (-oo, 0) . (5.1) 



Although the coordinates (3.4) cover only half of the fully analytically extended de Sitter spacetime, 
they exhibit the spatial homogeneity and isotropy and standard FRW form explicitly, and are convenient 
for admitting the standard Fourier mode decomposition in the flat spatial coordinate x. This allows 
the linear response equations to be expressed as ordinary differential equations in the time variable t 
(or r/), for each spatial Fourier mode separately. 
In the maximally symmetric de Sitter spacetime, 

R a c b d = H 2 (g ab g cd -5 a d 5 b c ) , (5.2a) 

R a b = 3H 2 5 a b , H = -, (5.2b) 

a 

R = 12H 2 , V a R = 0=aR, (5.2c) 

E ab = -2 (% ab = QH A g ab , E = 24H 4 , (5.2d) 

C a c b d = 0, A ab = B ab = 0, F = 0. (5.2e) 

Because of the 0(4, 1) maximal symmetry, with 10 Killing generators, the conformal "vacuum" of a 
conformal field in de Sitter spacetime also is maximally symmetric. The linear response approach 



requires a self-consistent solution of equations (1.1), around which we perturb the metric and stress 



tensor together. Hence the choice of this maximally symmetric state, known as the Bunch-Davies 
(BD) state is the natural one about which to consider perturbations. From the point of view of a free 
falling observer this conformal "vacuum" is in fact a thermally populated state with the temperature 
T H = H/2ir [3]. In this state 0(4, 1) de Sitter symmetry implies that (T a b ) is also proportional to 5 a b . 
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The self-consistent value of the scalar curvature R including the quantum contribution from (T a b ) is 
given by the trace equation, 



R + 4A = 8irG N b'E 



or 



H 1 



A 
3~ 



16irG K b'A 



+ 



(5.3a) 
(5.3b) 



in an expansion around the classical de Sitter solution with SttG^ A|6'| <C 1 [33 . Thus, in this limit 



the stress tensor source of the semi-classical Einstein's equations (1.1) in the BD state is a small finite 



correction to the classical cosmological term in the self-consistent de Sitter solution. The fractional size 
of this correction is hG N A./c' <C 1, if the de Sitter horizon scale H -1 is much greater than the Planck 



scale Lpi — 



1.6 1 6 x 10 33 cm. 



Given the smallness of the quantum correction in a self-consistent solution of (1.1), at first sight it 



would appear unlikely for the quantum fields to have any significant effects at the macroscopic scales 
of cosmology. However, these effects will become apparent only when both the geometry and quantum 
state of the matter are not fixed but allowed to vary dynamically. Since a quantum state specified 
over a macroscopic region of space is highly non-local, variations of the state are independent of the 
magnitude of local curvature variations. Although there is a very small coupling L 2 pi H 2 <C 1 controlling 
the loop expansion around the classical de Sitter background, the expansion in this coupling may be 
non-uniform or even infrared divergent due to state dependent variations on the horizon scale H . 
In the auxiliary field effective action approach, information about the quantum state of the underlying 



conformal matter fields is contained in the particular solution of the auxiliary field equations (4.6). In 



de Sitter spacetime the conformal differential operator A4 of (4.7) factorizes, 



A 4 |d 5 = -□(-□+ 2H 2 



(5.4) 



and the equation satisfied by ip is inhomogeneous whereas the equation for ip is homogeneous. Thus it 
is consistent to set ip = in order to obtain a de Sitter invariant state. Taking ip 7^ will not lead to a 
de Sitter invariant stress tensor. In any spatially homogeneous and isotropic state <p can be a function 
only of time. With ip = (p{t) one easily finds that the general homogeneous solution of either of Eqs. 



(4.6) is the linear combination, cq + c\a 1 + cia 2 + 03a 3 . All of these behaviors except the first lead 



to de Sitter non-invariant stress tensors, with the constant Co dropping out of the stress tensor entirely. 
Hence to obtain a de Sitter invariant state all the coefficients of these homogeneous solutions of Eqs. 



(4.6) must be set to zero. The remaining inhomogeneous solution to (4.6a) is easily found to be 2Ht. 
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If the solutions 



(p = 21nft = 2Ht 

$ = 0, 



are substituted into the stress tensors (4.4) and (4.5) the result is 



T%[(p, ^ = 0] dS = b'E a b [g dS ; <p] = 6b' H 4 5 a 



(5.5a) 
(5.5b) 

(5.6) 



which is exactly the value of the renormalized stress tensor of a conformally invariant field of any spin 



in the BD state, found either by direct computation or the methods of [26]. Thus the solutions (5.5) 
correspond to the Bunch-Davies state. 



With the self-consistent BD de Sitter solution (3.4), (5.1), (5.3), (5.5), and (5.6) one may consider 
the linear response variation, 

5R\ - ^6% = 8itG n 5 {{T a b ) loc + (r a fe )™} 

= SitG N { -a R 5A a b - f3 R 5B a b + b' 5E a b + b 5F a b } , (5.7) 

where all terms to linear order in 5g ab = h ab and in the variations of the auxiliary fields, 

Sep = ip — (p = 4> (5.8a) 
5^ = i;-i> = 'ip, (5.8b) 

are to be retained. All indices are raised and lowered by the background de Sitter metric g a b at linear 
order in the perturbations. 

For the first local tensor A a b it is simplest to vary the alternative form, 

A\ = -4C a cbd R cd -2UR\ + 2 -V a V b R + ^6 a b nR + 4R a c R c b - A -R\R - 5 a b R c d R d c + l -5\R 2 , (5.9) 
which follows from the Bianchi identities. The variations are simplified by using the symmetry properties 



of de Sitter space catalogued in (5.2), and by imposing the 5R = condition discussed in the previous 



section. It is easily demonstrated that this implies that 5(V a Vf,ii) = and 5(L\R) =0, as well, while 
5(DR a b ) = OSR a b . Hence, 



dS, SR=0 



□ + 



5R\. 



For the second local tensor B a b one finds, 



8B\ 



dS, 5R=0 



-2RSR a b . 



(5.10) 



(5.11) 
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The variation of the stress tensors obtained from the anomaly action, E a b and F a b depend on both 
the variations of the metric and the variations of the auxiliary fields. It is shown in Appendix B that 



5F\ 



dS, 6R=0 



2Ht 5AK 



IT. 

3 a? 



df + Hd t 



(5.12) 



From (4.15) one expects the variation of E a b to contain the local geometric term, 
- 25 ^H a h + \5B a h = -^ 5R\ - 2 



5R\ 



20H 2 



(5.13) 



3 9 3 
This has been verified by direct calculation using the algebraic manipulation programs Mathematica 



and MathTensor. The variation of the auxiliary field equations, (4.6) gives 



5{u\) - ^S(n<p) = (□ - f ) ti Uv) - -2(T,,vv) ( y/-r, 



*(nV)-?W) 

6 



(5.14a) 
(5.14b) 



Both this variation and that of the additional terms in 5E t t dependent upon 5ip = (ft can be simplified 
somewhat by a convenient choice of gauge. Consider 



$(□¥>) = 5(g ab V a Vw) = -h ab V a V b y - g ab (5r c ab )V c v +U<f>, 



,ab 



(5.15) 



where 



W C ab = 2 i-V c h ab + V a h c b + V b h c a ) 



g aD 6T 



ab 



Vah° 



-V c h, 



(5.16a) 
(5.16b) 



is the variation of the Christoffel symbol and its trace (with h = g ab h ab ). The fact that ip is a linear 
function of t allows the variation to be taken inside one wave operator, i.e. S(n 2 ip) = Ud(U<p) as in 



(5.14a), but from (5.15), 5(\Hf) ^ in general. However the additional terms in (5.15) may be set 



to zero by making the gauge choice, 



h ab V a V b (p = 



g bc 5Y a bc V a (p = I V b h a0 - ~V a h 1 V a <P = . 



ab 



1. 



(5.17a) 
(5.17b) 



Because of the simple form of (p from (5.5), these two gauge conditions are also equivalent to 



h ij gij = h a a + h tt = 

V*fcti = g ij ^ihj = \d t h tt 
2 



(5.18a) 
(5.18b) 



29 



or in terms of the decomposition (3.12), 



V 2 B = a{A + QHA) , 



C = 0. 



(5.19a) 
(5.19b) 



With this choice of gauge the variation in (5.15), 5(\3<p) = d<j>, and the auxiliary field Eqs. (5.14) can 
be written in Fourier space in the simple form, 



■ 5H '[ ■ (i// 2 • lr< ■ ! v - n. 
at z at 

^ + 5Hj t +6H 2 + k 2 e- 2Ui ) (<r-2h n ) . 



') 



(5.20a) 
(5.20b) 



with 



H 2 v 



H 2 w 



a 2 \drj 2 

1 (^ + k 2 
a 2 \ drj 2 



(5.21a) 
(5.21b) 



In terms of these quantities, with the condition 5R = and the gauge choice (5.17), the full variation 
of E\ is 



5E l 



20H 2 , 2 V 2 



d 2 + Hd t -—, 



3 a 2 



V' (w - 2h tt ) 



while that for F t in (5.12) becomes 



2H z h t 



8F\ 



2H 2 -> 

= AHt (- □ + AH 2 ) 8R\ =- V 2 v . 

dS,6R=o y 1 3a 2 



(5.22) 



(5.23) 



The gauge choice (5.17) is useful for simplifying the scalar auxiliary field contributions to (5.22) and 



(5.23), while the local geometric terms involving 5R t t in both expressions are independent of this gauge 



choice. In the next section it is shown that in fact the quantities w — 2hu and v are gauge invariant. 



All the variations in (5.10), (5.11), (5.22) and (5.23) involve the local geometric variation SR t . Thus 



the linear response equation in de Sitter space derived from the auxiliary field effective action takes the 
form, 



a 

H 2 



! (- □ + AH 2 ) 5R\ + (P B + 5e')5R\ + §(-□+ 4 # 2 ) SR \ 



+ ^ V 2 (w-2h tt ) 
2a z 



6a 2 



V v. 



(5.24) 
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where the coupling constants have been written in terms of the dimensionless small numbers, 

a R = 16nG N H 2 a R , (5.25a) 



192itG n H 2 P e 



327rG N H z b, 



Defining the dimensionless variable q by 

2a 2 



8G\ 



2a 2 
IP 



5R\ - -5R 



(5.25b) 
(5.25c) 
(5.25d) 

(5.26) 



it is straightforward to show that under the condition 5R = and with the gauge choice (|5.17|) 

v2 



(- □ + AH 2 ) 5R\ 



df + 7Hd t + 12H 2 



a" 



qH 2 

2a 2 " 



El 
2a 2 



df + 3Hd t + 2H 2 



(5.27) 



Thus, multiplying (5.24) through by —2a /H and Fourier transforming the spatial variables gives 



k 2 



with 



{1-0 R - 5e')q = eHtVq -a R Vq + e' j^(w - 2h tt ) - -j^v . 



(5.28) 



(5.29) 



Eqs. (5.20) and (5.28), with (5.21), (5.25), and (5.26) are the linear response equations for de Sitter 



space derived from the anomaly action and stress tensor. Since by assumption /? and e' are very small 



quantities they may be neglected with respect to unity on the left side of (5.28). Note that in the 



strict classical limit, where all the small quantities in (5.25) first order in H are set to zero, the classical 



constraint equation q = or 5R l t 



is recovered from (5.28). This is consistent with the absence of 



any dynamics in the pure Einstein theory in the sector of scalar perturbations without matter sources. 
This form can be compared now to those that are obtained by the method of Ref. jHl [10] reviewed 



9 



in the previous sections. From (1.3), (3.6) and (3.7), upon raising one index with the inverse metric 
ab - tt- 2 g ab , one finds 



SIR 1 



8irG K 



dr/K(rj - rf\ k; fx 2 ) SA^tf; k) 



'/ii 



(5.30) 
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with 



a R 5A\ - ( (3 R + - ) 8B\ - 2b' 5 {3) H 



ft 4 



ti j 
t 



(5.31) 



Upon making use of fl418| ), (|oTT0|), (|5T27|), 0^26], p^ty, (|5T25]) and §iJty, Eq. <|5T30]) becomes 

5e^ 



21^ 
IP 



eHt Vq-a R Vq- ^ J dr( K( V - rj- k; M ) [& Pg]^ - ^<5(T' 



(5.32) 



Comparing to (5.28), it is apparent that the linear response equation (5.32) based on the methods of 



[21 [TU] differs from that based on the anomaly action in four respects: 



The coefficient of SB t t in (5.30) and (5.31), — (5 R — ^ is r eplaced by -0 R + | in M and p.28\. 



This is an inessential difference, amounting to a different coefficient of the finite local R 2 term in 



the action, already discussed, and which, as can be seen from (5.25), is any case negligible. 



The linear response form (5.30) or (5.32) contains the general state dependent term, 5(T a b ) R , 



whereas (5.28) contains the specific w — 2hu and v auxiliary field terms. 



The auxiliary field action gives rise to additional equations of motion (5.20 ) for the state dependent 
terms, absent in the approach of [HI ITU] . 



Linear response based on the anomaly action, (5.28) lacks the term with the non-local kernel K 



in (5.30) and (5.32). 



In O [10] , the variation of the quantum state of the field was (implicitly) constrained by the metric 
variation, and not allowed to vary independently. This is equivalent to setting v = w — 2hu = in 



the anomaly auxiliary field approach. Since these fields obey the homogeneous equations (5.20), the 



particular solution in which they both vanish is allowed. However, the local auxiliary field approach 
allows for a wider class of variations of the state than that of [9j [TO] . That the w — 2htt and v auxiliary 



field terms in (5.7) are the result of allowing the state to vary over and above the local metric variation 



follows from the fact that the freedom to vary the auxiliary fields amounts to freedom to vary the 
boundary conditions on the Green's function of A4 entering the non-local form of the trace anomaly 
effective action |16j . From (5.24) one sees that these terms scale with the conformal factor exactly as 



expected for the state dependent variation S(T t t ) R in the general form (5.32). In (5.7) these terms take 
a specific form related to the variation of the scalar auxiliary fields which obey their own independent 
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equations of motion (5.20b) and (5.20a). The important feature of the auxiliary field formulation of 



the anomaly effective action (4.1)-(4.2) is that these additional state dependent terms are expressed in 



terms of local field variations that have their own independent equations of motion (5.14). Although 



such state dependent variations are certainly allowed on general grounds, the specific form of these 
variations are determined by the local auxiliary field form of the effective action. 

VI. GAUGE INVARIANT VARIABLES AND ACTION 



As already remarked in Sec. Ill the scalar sector perturbations involve the four functions A, B,C, £, 
only two linear combinations of which are gauge invariant. In this section gauge transformations are 
discussed and a set of gauge invariant variables are constructed. The linear response equations around 
de Sitter space are written in terms of these gauge invariant variables, and the gauge invariant quadratic 
action functional corresponding to these equations is given. 

The linearized coordinate (gauge) transformation of the metric perturbation is 



h a b — ► h a b + V a ^fe + Vb£ a . 



Under this gauge transformation the scalar auxiliary fields transform as 



(6.1) 



(6.2a) 
(6.2b) 



With the definitions 



e = aT, 

C = 7] lj djL . 



(6.3a) 
(6.3b) 



appropriate for scalar perturbations, it is found that the linearized metric functions in the decomposition 



(3.12) transform as 



A^A + tiT + aT, 

B —> B + aL - T , 
1 ->•> 

C + - V L + aT . 
£ -^£ + L . 



(6.4a) 
(6.4b) 
(6.4c) 
(6.4d) 
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These are the linearized coordinate transformations possible in the scalar sector |29l [30] . By considering 



the first order variations ofEqs. (6.2) and making use of the notation introduced in Eq. (5.8), one obtains 



the transformation for both auxiliary field perturbations: 

<j)^(j) + 2aT , 
ip — ► ifi . 



(6.5a) 
(6.5b) 



Thus ip is already gauge invariant, as is the quantity v defined by (5.21a), while <\> transforms due to 



the non-zero tp = 2Ht in the BD state. It is easily checked that the from the transformation of the 



field (6.5a) that 



■ + 2aB - 2aa£ 



(6.6) 



is gauge invariant. This is similar to the gauge invariant variable that can be constructed from the 
scalar field in scalar inflaton models of slow roll inflation [34 . Finally the metric variables, 



T A = A + dt(aB)-dt(a 2 d t £), 



V 2 £ 



+ aB — ad£ 



(6.7a) 
(6.7b) 



are invariant under the linearized gauge transformations (6.4). These are the gauge invariant Bardeen 



Stewart potentials denoted by $ A and <J?c i n Refs. [2§] EUJ 

The two quantities 5R and q encountered in the linear response analysis can be written in terms of 
the metric gauge invariant variables T A and T^. Indeed, using the variation for 5R [35] one finds 



5R 



■nh + V a V b h ab - R ab h 



ab 



6(T C - HT A ) + 2AH{T C - HT A ) 



V 2 T C , 



(6.8) 



Hence condition (3.11 ), 5R = 0, is gauge invariant, and provides one constraint between the two gauge 



invariant potentials T A and T^. Next one can verify that q, defined in Equation (5.26), can be written 
in the form, 



2a 



fp 5G\ = 12a<l-r c -r A 



1 



H 2 



V 2 T C , 



(6.9) 



which is also gauge invariant. 

In the linear response equations the quantity w — 2hu is also encountered in the gauge defined by 



the conditions (5.19). The gauge invariant quantity which reduces to w — 2hu in this gauge is 

~d 2 



H z u 



dt 2 + dt a 2 



<£> + 6H 



d ^-2H d ^-8H 2 T A . 
dt dt A 



(6.10) 
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This can be seen by evaluating Ya, Tq, and $ in the gauge (5.19) using the definitions (6.6) and (6.7). 

Then the linear response equations around de Sitter space can be written in terms of the gauge 
invariant variables u, v, SR, and q. With the condition 5R = imposed, the equations for u and v are 

„2- 



2 
dt 2 



dt 
d 



+ 5H-+6H Z + 
dt 



0, 
0, 



(6.11a) 
(6.11b) 



while that for q is 



2& 

= eHt Vq — a R Vq 



^ dr } 'K{r ] -r 1 ';k;ii)[n 2 Vq) nl + e 



,k 2 

w 



(6.12) 



Here the differential operator D is defined by (5.29). In (6.12) the (3 R , e and s' terms on the left 



side of (5.28) or (5.32) which differ in the two linear response approaches but which in any case are 
small compared with unity, have been dropped. The non-local term from the analysis of [§1 [TU] has 
been retained and the specific state dependent terms from the auxiliary field anomaly action have been 



written in terms of the gauge invariant variables u and v. In arriving at (6.12) we have reconciled 
the exact but formal calculations of [9j [10] with those following from the anomaly effective action by 
adding the non-local term from the former to the specific state dependent variations determined by the 
auxiliary fields of the latter approach. By the arguments of Sec. Ill, all the gauge invariant information 
about scalar perturbations around the self-consistent de Sitter invariant BD state and all components 



of the curvature tensor variations can be obtained from the solution(s) of (6.11 )-(6.12 ). This is the final 
gauge invariant form of the linear response equations for the only remaining metric degree of freedom 
in the tracefree but spatial scalar sector of perturbations of de Sitter spacetime whose solutions we will 
analyze in the next section. 

If the classical Einstein-Hilbert action is added to the effective action of the anomaly, Eqs. (4.1) and 

(4.2) , and both are expanded to quadratic order about the self-consistent BD de Sitter solution given by 

(5.3) , (5.5), and (5.6), the resulting quadratic action can be expressed in terms of the gauge invariant 
variables u, v, T^, and Tc in the relatively simple form, 



[1 - be') S G + b' J d 3 xdta 3 j- 



H 4 u 2 H 2 u5R 

— — + 



d 3 x dt a 3 



H 2 v5R 41na 
+ 



3 a 4 



V 2 (T A - T c ) 



(6. 
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where 



S G 



8ttGn 



d 3 x dt a 6 



dT 



c 



Of 



+ QHT A & ^ + \ (VT A )-(VT C ) + 
at a z 



(VT C ) 2 



3H 2 T 2 A 



(6.14) 



is the Einstein-Hilbert part of the action, and 5R is given by (6.8). Varying (6.13) with respect to 
<&,iP,Ya and Tq, and setting 5R = yields the gauge invariant linear response equations (6.11) and 
(6.12), without the non-local contribution involving K, and without the a R and (5 R terms which would 
require adding to (6.13) the contributions of the purely local C 2 and R 2 actions also expanded to 
quadratic order in the perturbations about de Sitter space. 



VII. SOLUTIONS TO THE LINEAR RESPONSE EQUATIONS 
A. Solutions of the first kind: u = v = 



Since u and v satisfy the same homogeneous equation, cf. (6.11 ), we first consider the case u = v = 0. 



Eq. (6.12) is then homogeneous in q and includes only those variations of (T a b ) R which are driven by 
the metric fluctuations, i.e. variations of the first kind. The homogeneous equation possesses the trivial 
solution q = 0. This is the unique classical solution in pure classical gravity, for when h = 0, all effects 



of the stress tensor (T a b ) R and its quantum fluctuations are set to zero on the right side of (6.12), and 



there is no dynamics at all in the scalar sector. The classical equation, q = is an equation of constraint. 



Next it is straightforward to find non-trivial solutions of (6.12) if u = v = and the non-local term 
involving K is neglected. One may then check the consistency of this local approximation by inserting 
the previous solution into the integral and evaluating the non-local term. Since at late times the eHt T>q 
term dominates over the 6lr T>q term, which in any case just shifts the origin of the time variable t, the 
cur T>q term can be neglected as well without loss of generality. 



If the non-local term in Eq. (6.12) is ignored and u = v = = aR = then the linear response 
equation reduces to 

q = eHtVq (7.1) 

The Fourier components with k ^ are red-shifted away exponentially rapidly by the de Sitter expan- 
sion, k 2 e~ 2Ht — > 0. Hence after a few expansion times this term may be neglected, effectively setting 
k = 0. Then, 
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This equation describes an exponential growth on the time scale H/y/e ~ Mpi. Substituting the 
exponential form, 

cHt 



q(t) = exp (^J 



W{t)<1t 



(7.3) 



one finds 



1 dW 



+ W + 3W + 2 



1 



H dt eHt 
Ignoring the dW/dt term, the zeroth order WKB solution is 



(7.4) 



W 



3 1/ 4 

-2 + 2i l + jm 



(7.5) 



Starting from an initial time to such that Hto is of order unity, the WKB approximation (7.5) is valid 
up until the turning point where W vanishes. This turning point is at Ht m ax = l/2e and this is where 
q achieves the maximum value, 

'In 2 2Ht ' 



i, 



£ < 1. 



(7.6) 



£ y/£ 

Here to is the initial time which is assumed to be much less than (i?e) _1 . Thus if the non-local term 
is ignored, then the perturbation becomes exponentially large in e - 1 ~ MlJH 2 » 1. The additional 
factor of the time t reduces this growth at later times, and after a time of order e _1 eventually turns 
off the exponential growth. 



Eq. (7.1 ) was also solved numerically for e = 0.01 and the results are shown in Figure ??. Note that 



the solutions grow rapidly for a long time before reaching a maximum and finally decreasing again. The 



exponential growth is well described by the analytic WKB approximation of (7.3) and (7.5) 



To test whether this general behavior survives in the full linear response equation including the 



non-local term, one can first solve (7.1 ) and then substitute the solution into the non-local term to see 



if its effects are significant or not. The non-local term in Eq. (6.12) can be written in the form, 

r-ij 



1 



2ft 2 



drj'K(r] - rf; k\n)f{rf) 



(7.7) 



with K given by (3.5) and 



f( V ) = en 2 (r))Vq(ri) 



(7.8) 



The behavior of the kernel K is discussed in Appendix A. For an initial value formulation of linear 



response, the lower limit of the time integral in (7.7) should be taken as an arbitrary but finite 7/0 . 
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FIG. 2: The plot shows solutions to the linear response equation (6. 12 1 in the approximation (7.1) when the 



non-local term is neglected, u — v — 0, e — 0.01, and k — 100. The initial time is Ht = 10. The solid curve is 
the solution which begins with q = 1, q = and the dashed curve is the solution which begins with q — 0, q = 1. 



The solutions grow exponentially according to (7.3) and (7.5 1 up until Ht max ~ ^ = 50. 



As discussed in Appendix A, a proper regularization of the non-local kernel K, such as by means of a 
Pauli-Villars regulator mass, produces transient terms which fall off at late times, and can be neglected. 



Then Eq. (A13) can be used with the lower limit of the integral replaced by rjo -C rj to obtain the 



approximate form, 



no 



(7.9) 



For the solution of (7.1) 



n 2 (r])q(r]) 



- In(-Hri) 



exp 



W(t)oIt 



(7.10) 



ln(n(r?)) H 2 r] 2 ln(-Hr]) 

In Figs. [3] and [4] the non-local and local terms on the right hand side of the linear response equations 
are shown for the case in which q is the solution to the local linear response equations displayed in Fig. 
[2] It is clear from the Figures that for the solutions shown the contribution of the local and non-local 



terms on the right hand side of (6.12) are not only comparable but tend to partially cancel. Hence we 



conclude that it is not correct to ignore the non-local term relative to the local terms in (7.1) for these 
homogeneous solutions for which u = v = 0. 

It is possible to understand the behavior of the non-local term analytically for the case that f(rf) is 



a rapidly increasing function, as it is for (7.10). In this case the largest contribution comes from the 



region near the endpoint of the integral (7.7) as rf — > r]. Consider a time rji < r\ close to the upper limit 
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FIG. 3: Shown are the non-local and local terms for the case k = 100. The dashed line corresponds to the local 
term eHtDq — q. The initial conditions for the solution shown for q are specified at Ht = 10. 
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FIG. 4: Shown is the approximate cancelation which occurs when the local and non-local terms, eHtDq = q and 
I respectively are added together for the case k = 100. 

of the non-local integral, so that the corresponding interval Ai] = rj — rji satisfies 



kAr] < 1 . 



(7.11) 



If this condition is satisfied then from (All ), c\[k{r] — r/)] ~ C + ln[/c(r/ — rf)] and (7.9) becomes 



drj' In (fe( 7? - r /))^M + 



Cf(r,)-Cf(n-Ar,) + f(ri)\n 



^ {[/(»?) " f(v ~ [HkAr,) + C]+ f(r,) In 

{ mn) Al]d ^r [HhA7]) + c] + f{7]) ln 111 



(7.12) 
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Now from (7.3), (7.5) and (7.10), during the period of exponential growth, 

/ 



Hence if 



d>V ~ \v\ ~ \ri\y/e\hi(-Hrf)\ 



Ar] < \r]\^e\ln{-Hr])\ , 



(7.13) 



(7.14) 



for which (|7.11|) is well satisfied, with e<1, then 



h 



< eVqln ( i—e 



-e{Ht)Vq + eVq In ^e c V£ffi) 



(7.15) 



The first term of (7.15) exactly cancels the local term, e{Ht)T>q on the right hand side of the linear 



response equation (6.12), thus giving an analytic approximation which reproduces the cancelation found 



numerically, and exhibited in Fig. [4] The second term in (7.15) gives an estimate for the remainder. 
Examination of the steps leading to ( |7.15 ) shows that this approximation is valid as long as T>q is a 



rapidly increasing function. In this case, the non-local integral (7.7) receives its dominant contribution 



from the very short time ultraviolet region (7.14) close to its upper endpoint. 



From both this analytic estimate and the numerical studies we conclude that for the homogeneous 



solutions of the linear response equation (6.12), the non-local term cannot be ignored. Thus the expo- 



nential growth found in (7.3)-(7.6) by neglecting the non-local integral in (6.12) is not reliable. Actually 



this could have been anticipated by returning to the origin of these terms in Section III. The non- 
local term, and the local a R T>q and HtDq terms, all of which involve T>q have the same origin in the 



logarithmic distribution in Fourier space. Eq. (2.8) shows that the T>q dependent terms should be 



considered together, since only the sum is independent of the arbitrary renormalization scale \x. The 
Ht Dq = In 0, T>q term arises from the fact that the logarithmic distribution is defined in flat conformal 
coordinates, in terms of the conformal frequency and momentum K a = (oj,k), whereas the frequency 



and momentum relative to the physical line element (3.4) is K a /Q (for slowly varying O). Thus one 
should expect the combination, 



In 



-uj 2 + k 2 — iesgnu 
V 2 



21nft = In 



-uj 2 + k 2 — iesgnuj 

^2 



(7.16) 



always to appear together. Heuristically, the local In £1 T>q term is just the term needed to insert fl in 



the denominator of the logarithm in (7.16), and convert the conformal frequency and wave-number to 
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their values in the physical metric. Since it is known from the linear response equations in fiat space 
that the only non-trivial modes occur at the Planck scale, these same Planck scale solutions should 



persist in the de Sitter background when all three of the terms involving T>q in (6.12) are considered 
together and solved self-consistently. For these homogeneous solutions of (6.12) for which u = v = 



the non-local integral / cannot be neglected. This is the significance of the integral being dominated by 
its extreme short distance logarithm near its upper limit. To find these modes accurately one should 



solve the full non-local Eq. (6.12). However, these ultraviolet Planck scale solutions should be excluded 



from consideration in de Sitter space in any case for the same reason as in flat space, because they lie 
outside of the range of validity of the low energy semi-classical description of gravity. Therefore we do 



not pursue the solutions of the homogeneous Eq. (6.12 ) further, and turn instead to the inhomogeneous 
solutions of the second kind. 



B. Solutions of the second kind: u ^ or u ^ 



Next consider the case in which the inhomogeneous gauge invariant state dependent terms u or v 
are different from zero. The general solution of (6.11) for either u or v is easily found. In either case 
the solutions are 



1 



u± = v± 



, exp I ±-^- ] e tk ' x = rj 2 exp(=Fi/c?? + ik ■ x ) 
H z a z \ Ha J 



(7.17) 



Note that this solution and the Eqs. (6.11) it satisfies involve the cosmological horizon scale H, but not 
the Planck scale. Thus we term these new solutions cosmological horizon modes. Taking e.g. u = 0, 



neglecting all the T>q terms in (6.12), we have 



e -9 



q ~ -p H 2 r) 2 e Tikr i +ik -* 



This corresponds to a linearized stress tensor perturbation of 

T2 „ k 2 H 2 



5(T l 



H a 

" ^ ^ 11 ^fikri+ik-x 



16irG N a 2 



(7.18) 



(7.19) 



The solutions with u of the form (7.17) give linearized stress tensor perturbations similar to (7.19) with 
b replaced by b' . 



We will now show that it is legitimate to neglect the terms involving T>q in (6.12) for these solutions. 
Note that 



Vq = ±— q = ±2i^~ q ■ 
Ha H 



(7.20) 
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where k p h ys = k/a = H\krj\ is the redshifting physical momentum of the mode. Recalling the definitions 



(5.25 ) we have 



provided 



\a R Vq\ ~ \a R G N Hk phys \ \q\ <C \q\ 



G N H 2 ^1, or H^Mpi, 
G N kp hys «1, or k phys < Mpi 



(7.21) 

(7.22a) 
(7.22b) 



and \a R \ is of order unity. Thus, from (7.21) the a R T>q term in (6.12) can be neglected. 



For the non-local term note that (7.20) falls off with time, which has the consequence that the 



function (7.8) appearing in the integral also falls with time. In this case the nonlocal integral / in 



Eq. (7.7) receives contributions over the entire cosmological expansion time scale H 1 , and remains 



bounded. This is in contrast to the non-local term evaluated on the steeply rising function considered 
in the homogeneous case when u = v = 0, where (7.7) was dominated by the short time behavior (7.14) 



near its upper endpoint. For the solution in (7.18), / may even be calculated analytically using the 



Pauli-Villars form of the kernel in Eq.(A21), for M large and setting M = fi in the terms that remain 



after this limit is taken. In this way one finds that all the terms in the non-local integral are small and 



may be neglected, provided the conditions (7.22) are satisfied, except for one possibly large logarithmic 
term which is of order, 

,2 



bG^H- In 



l 1 
Hk 



Because of (7.16) with $7 = a, this combines with the remaining local term involving T>q, i.e. 

k 



e HtVq ~ bG N H- ln(a) q 



to give 



bG„Hh 



phys 



In 



(7.23) 



(7.24) 



(7.25) 



H kphy S , 

which is also of negligible magnitude with respect to l^l provided k p hy S < ^ < Mpi, b is of order unity, 



and conditions (7.22) are satisfied. Since these are necessary conditions for the applicability of the 



semi-classical effective theory (1.1) in the first place, we conclude that all the T>q terms in the linear 



response equation (6.12) may be neglected compared to (7.18), and (7.18) is a non- vanishing solution to 



the full linear response equations, to a very high degree of accuracy if the conditions (7.22 ) are satisfied. 
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The other components of the stress tensor perturbation can be found for this solution from (7.19) by 



a general tensor decomposition for scalar perturbations analogous to Eqs. (3.12) for the metric. That 
is, the general perturbation of the stress tensor 5{T ab ) R in the scalar sector can be expressed in terms 
of 8(T U ) R plus three additional functions. These three functions are determined in terms of 5(T a ) R by 
the conditions of covariant conservation, 

V b 5(T ab ) R =0 (7.26) 
for a = t and a = i (two conditions), plus the tracefree condition, 

S(T a a ) R = 0, (7.27) 
imposed as a result of the 5R = condition. A straightforward calculation using these conditions and 



the Christoffel coefficients (3.16) gives then 



k 2 - hH 2 -* 

5(T a ) R = bH 2 -^ e ^kv+ik-x = -— V 2 e Tikv+ik-x ) ( 7 _ 28a ) 

Uu - hH 2 f> 2 

5(T U ) R = ±bH 2 ™ e Tikr,+ik-x = v_n_ c, eTikv+ik . x ^ (7 2gb) 

hiuj - hH 2 f) 2 

5{T ij ) R = bH 2 e* ikr > +ik - x = —j— e ^ +ik - x . (7.28c) 

CL CL OX OXJ 

for the other components of the stress tensor variation for these modes in the flat FRW coordinates of 
de Sitter space. 



Thus, the auxiliary fields of the anomaly action yield the non-trivial gauge invariant solutions (7.28) 
for the stress tensor and corresponding linearized Ricci tensor perturbations 5R ab . Being solutions of 
(6.11) which itself is independent of the Planck scale, these solutions vary instead on arbitrary scales 
determined by the wavevector k, and are therefore genuine low energy modes of the semi-classical 
effective theory. The Newtonian gravitational constant G N and the Planck scale enter Eq. (6.12) only 
through the small coupling parameters e and e' between the auxiliary fields and the metric perturbation 
q. Thus in the limit of either flat space, or arbitrarily weak coupling G N H 2 — > these modes decouple 
from the metric perturbations at linear order. 



The result (7.28) would be obtained if the anomaly action (4.1 ) alone is used to generate the linear 
response Eqs. (5.28), and the non-local term is neglected completely. This demonstrates the relevance 
of the anomaly action for describing physical infrared fluctuations in the effective semi-classical theory 
of gravity, on macroscopic or cosmological scales unrelated to the Planck scale. 
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VIII. COSMOLOGICAL HORIZON MODES 



The modes found in Sec VII B when the gauge invariant variables u and/or v are nonzero can vary 



on any scale rather than the Planck scale, which characterizes modes of the first kind found in Sec VII A 



when u = v = 0. The second set of modes arise naturally from the scalar auxiliary fields which render 
the non-local anomaly action local, and are therefore implied by the form of the trace anomaly at 
one-loop order. This is a non-trivial result since these modes appear in the tracefree sector of the 
semi-classical Einstein equations, with 5R = 0, and hence cannot be deduced from the form of the trace 
anomaly itself, but only with the help of the covariant action functional and the additional scalar degrees 
of freedom which the anomaly implies. These additional modes, which couple to the scalar sector of 
metric perturbations in a gauge invariant way, are due to a quantum effect because the auxiliary scalar 
fields from which they arise are part of the one-loop effective action for conformally invariant quantum 
fields, rather than a classical action for an inflaton field usually considered in inflationary models [34 . 

In Section V it was shown that the background solutions (p = 2 In $7 = 2Ht and ip = result in the 
stress-energy tensor for conformally invariant fields in the Bunch-Davies state. The additional modes 
arising from perturbations of the auxiliary fields tp and tp from their background values correspond to 
changes of state for the underlying conformal quantum fields from their de Sitter invariant BD state. 
Some physical intuition about these modes and the changes of state of the underlying quantum fields 



they correspond to may be gleaned from the form of their stress tensor in (7.28). If one averages this 
form over the spatial direction of k, a spatially homogeneous, isotropic stress tensor is obtained with 
pressure p = p/3. Thus in FRW coordinates this averaging describes incoherent or mixed state thermal 
perturbations of the stress tensor which are just those of massless radiation, which redshift with a -4 . 

A second interpretation of this second set of modes emerges if we consider coherent linear superpo- 
sitions of different k solutions in different global coordinate systems. In Ref. [16] a class of solutions 



to the auxiliary field equations (4.6) were found for de Sitter space in static rather than cosmological 
coordinates. In that background field calculations the spacetime geometry was held fixed and there 
is no restriction on the state other than that required for renormalization of the stress tensor for the 



quantum fields. In the linear response equations the solutions to the perturbed equations (6.11) for 



the auxiliary fields become source terms for the metric fluctuation q in Eq. (6.12) so their backreac- 
tion effects on the spacetime geometry are included to linear order in fluctuations about the de Sitter 



background geometry. However, in either case the equations (4.6) for the gauge invariant perturbations 
of the auxiliary fields, u and v, being linear, are the same equations whether or not the background is 



44 



varied. Thus they correspond to the same changes of state as in the background field calculations of 
|16j , and can be compared directly with the fixed de Sitter background solutions found for the auxiliary 
fields in the static coordinates in |16| . 

To make this comparison we first need to consider the relationship between the flat FRW and static 
coordinate systems in de Sitter space. The de Sitter line element in its static form, 

dr 2 



ds z 



-(1 - H 2 r 2 )dt z + 



1 - H 2 r 2 



+ r 2 (d6 2 + sm J 



i.l) 



is identical to (3.4 ) with the de Sitter scale factor (5.1 ) if the static time t and radius vector r are related 



to the flat FRW coordinates (t, x) of (3.4) by 



\x\ e m 



pe 



Ht 



i-^ln(l-tfVe^) 



3.2a) 
i.2b) 



The inverse transformations are 



r e 



-Ht 



P = \ x \ 



t = t + 



2H 



Vl - H 2 r 2 ' 
ln(l -H 2 r 2 ). 



3.3a) 
3.3b) 



The Jacobian matrix of this 2x2 transformation is 



/ di 


at \ 




at 


dp 




I dr 


dr 




\ at 


dp ) 





1 Hr 2 
1-H 2 r 2 p{l-H 2 r 2 ) 

Hr r - 



(8.4) 



Using these relations, one may express the action of the differential operators in Eq. (6.11) in terms of 



the static coordinates (8.1) instead. We consider the case that the scalar auxiliary fields are functions 



of r only and focus attention on v. A short calculation using (|8.4|) shows that 

1 d 



V 2 . 

-f v(r) 



> dv\ 



2 dr \ dr J 



while 



zj2 / n ( d 2 , „d V|\ u2 2 1 d ( 2 d^{r) 

H v{r) =\W + H 8i-^ ^ = " (1 " H r ) V^r \ ~dr 



(8.5) 



(8.6) 



operating on functions v = v(r) which are functions only of r and not i. 
In [TBI the static solutions to 



A 4 V(r) = 
45 



(8.7) 



are given. The four solutions are 



1 



In 



1 - Hr 



In 



1 - Hr 



1. 



1 



Hr V 1 + Hr J \ 1 + -HV / <" 

The last of these is singular at the origin and so was not considered in [16]. In any case this solution 



and the third constant solution to (8.7) give vanishing contribution to v in (8.6) while from (8.6) the 
first and second solutions give for v, 



4 



4 1 + H 2 r 2 
Hr 1 - H 2 r 2 ' 



(8.9) 



respectively. The second gives a singular contribution to v and the stress tensor at r = 0, and we 
consider only the first solution, taking 

1 (8.10) 



1 - H 2 r 2 



Then, from fl5.23|), fl5.24|), (|8.5|) and fl8.10|), we obtain 

4H 4 



6(T l 



2bH 2 Vf 



+ 



16H 6 r 2 



8R\ = 8ttG N 5(T t t ) R = -sH 2 j 



(1-H 2 r 2 ) 2 3(1 -H 2 r 2 ) 3 
1 4H 2 r 2 



+ 



111a) 
111b) 



(1-H 2 r 2 ) 2 3(1 -H 2 r 2 f 

The solutions for u are exactly analogous. This shows that a linear superposition of solutions of the linear 
response equations of the second kind in static coordinates can lead to gauge invariant perturbations 
which diverge on the de Sitter horizon. 
To see what 



11) corresponds to in the static (t,r) coordinates, we use the form of the other 



components in (7.28) in FRW coordinates, 



8{T U ) 



and the transformation relation for tensors 

2 



bH 2 

,2 







+ 2H 



bH 2 d 2 v 

a 2 dx 1 dxi 



dv 



T 



tt 



with ( |8.4| ), ( |8.5[ ), to obtain 

= bH 4 



^ dt) ^\dx l 



-(I- H 2 r 2 )5(T l 
1 d f r dv ) + 



rptl _|_ 



112a) 
112b) 

(8.13) 



r dr \ dr 
16 b H 4 
(l-H 2 r 2 ) 2 



bH 4 



1-H 2 r 2 ) 



1 dv 
r dr 



5H 2 r 



dv 
dr 



.14) 
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Here use has been made of the following identities, 

dt dt dp dt Xi 



dx l 




dp dx l dp p 



dr d dt d 

at = ~dldr~ + dtdl 



d 



Xid 1 = x' l di = p 



d_ 
dp 



= Hr 

dr 

I dr d di d 
f \dpdr + d~pd~t 
d 



dr 



XiXid l d J = (x l di)(x J dj) — x l di = r— | r , 

dr \ dr 



dr 



valid when operating on functions of r only. Likewise we find 



6{T 



r \ 
VI R 



5{T e t 



5{T 4 



lQbH 4 



5.15a) 
5.15b) 
3.15c) 
L15d) 

(8.16) 



3(1 -H 2 r 2 ) 2 ' 

corresponding to a perfect fluid with p = p/3, in static coordinates. The contributions proportional to 
b' instead of b from the u solutions are also of the same form. 



The form of the stress tensor (8.14), (8.16) is the form of a finite temperature fluctuation away 



from the Hawking-de Sitter temperature T H = H/2tt of the Bunch-Davies state in static coordinates 



3d] . Since the equation the solutions (8.8) satisfy is the same as (4.6), it follows that there exist linear 



combinations of the solutions (7.17) found in Sec. VII B which give (8.10) and the diverging behavior of 



the linearized stress tensor on the horizon, corresponding to this global fluctuation in temperature over 
the volume enclosed by the de Sitter horizon. Note that in static coordinates the stress tensor p = p/3 
does not require averaging over directions of k, but a particular coherent linear superposition over modes 



(7.17) with different k in order to obtain a particular isotropic but spatially inhomogeneous solution 



of (5.20a), which selects a preferred origin and corresponding horizon in static coordinates (8.1). The 



fluctuations in Hawking-de Sitter temperature thus preserve an 0(3) subgroup of the de Sitter isometry 
group 0(4, 1). 



To follow the diverging behavior (8.14), (8.16) all the way to the horizon one would clearly require 



a linear combination of the solutions (7.17) with large Fourier components. However once 8irG N times 



the perturbed stress tensor in (8.14) becomes of the same order as the classical background Ricci tensor 
H 2 , the linear theory breaks down and non-linear backreaction effects must be taken into account. This 
occurs at r = H -1 — Ar near the horizon, where 



Ar ~ Lpi , 



or because of the line element (8.1 ) at the proper distance from the horizon of 



H 



(8.17) 



.18) 
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Thus in the background field calculation, one finds that the stress-energy is relatively small well inside 
the horizon. This correspond to a maximum k p hy S ~ 1/1 <C Mpi, where the semi-classical description 



may still be trusted. At the distance (8.18) from the horizon, the state dependent contribution to the 
stress-energy tensor becomes comparable to the classical de Sitter background curvature, the linear 
approximation breaks down, and non-linear backreaction effects may be expected. This is the same 
estimate that was obtained in Refs. [37], where a matching of de Sitter space within the horizon onto 
an exterior Schwarzschild spacetime was discussed. 

IX. SUMMARY AND CONCLUSIONS 

In this paper the specific form of the linear response equations in a de Sitter space background for 
semi-classical gravity coupled to conformally invariant matter fields has been derived with the principal 



result being Eqs. (6.11 )-(6.12 ). In terms of a decomposition of the perturbations into scalar, vector, 



and tensor parts with regard to spatial hypersurfaces |29[ 130] , the analysis has been restricted to the 



scalar sector in (3.12), and it has been shown explicitly that the resulting linear response equations 



can be written in terms of gauge invariant variables through (6.7)-(6.10 ). The quadratic action in 
terms of these gauge invariant variables corresponding to the linear response equations is given by Eqs. 
(6.13) and (6.14). A general condition, 5R = 0, has been utilized which eliminates conformal or trace 
solutions which vary on the Planck scale and are thus outside the range of validity of the semi-classical 
approximation. This condition is approximate in generic spacetimes but is an exact restriction on 
solutions in the maximally symmetric de Sitter spacetime. Remaining under consideration are only 



those solutions of the linear response Eqs. (1.3) around de Sitter space which both have zero four 
dimensional trace, but which are scalars with respect to the FRW spatial sections. Owing to the 
freedom to make linearized coordinate transformations of the background de Sitter space, this leads to 
one and only one remaining gauge invariant scalar metric degree of freedom, whose dynamics is fully 
described by the tt component of the perturbed semi-classical Einstein equations in FRW coordinates 



(3.4|), namely ( |6.12[ ). 

The importance of the trace anomaly and the auxiliary fields used to render it local have been 
emphasized in our analysis. Among other things, the linear response analysis presented here provides 
an important test for the auxiliary field action. If the auxiliary field form of the anomaly effective 
action is used, then the linear response equations consist entirely of local linear partial differential 
equations with state dependent perturbations resulting from variations in the auxiliary fields. In the 
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linear response equations derived from the exact effective action there are minor differences associated 
with the values of renormalization parameters, along with two more significant differences. One is that 
in the anomaly action approach the unspecified state dependent terms allowed in the general approach 
take quite specific forms in terms of the scalar auxiliary fields, which possess their own dynamical 
equations of motion. The second significant difference is that the anomaly action takes no account of 
the conformally invariant term associated with the non-local term with the kernel K. 

The solutions to the linear response equations have been investigated in Sec. VII. Our main results 
can be divided in two parts, depending upon whether the perturbations of the stress tensor are driven 
by changes in the metric (solutions of the first kind) , or additional state dependent perturbations such 
as those generated by the auxiliary fields of the anomaly action (solutions of the second kind) are 
considered. For perturbations of the first kind, the local linear response equation has been solved by 
first ignoring the non-local term involving K and then evaluating the non-local term using the solutions 
so obtained. It was found that the non-local term cannot be neglected for this first class of solutions to 
the linear response equations, and this is closely related to the fact that these modes are fundamentally 
short distance Planck scale modes, analogous to those already found in flat spacetime, which in any 
case lie outside the range of validity of the semi-classical effective theory. 

The only remaining solution to the exact linear response equations for perturbations of the first 
kind is the constrained trivial solution q = 0. Since only this trivial solution to the linear response 
equation survives when the Planck scale solutions are eliminated, it means that our criterion for the 
validity of the semi-classical approximation in de Sitter space is satisfied [6 j . As originally formulated 
this validity criterion applies only to gauge invariant fluctuations of the first kind, which depend upon 
the two-point retarded correlation function for the stress tensor, ([T a b , T cd ]) and thus probe the response 
of the geometry to the quantum fluctuations of the stress-energy tensor about its mean value. 

For state dependent perturbations of the second kind the situation is much more interesting. A new 



class of infrared modes has been found in the tracefree but scalar sector, given by Eqs. ( 7.17 )-( 7.19 ) 



and (7.28), which are associated with the scalar auxiliary fields of the anomaly effective action. Since 
these modes vary not on the Planck scale but on arbitrarily large scales which can be comparable to 
the cosmological horizon scale in de Sitter space, we have termed them cosmological horizon modes. We 
emphasize that the existence of these scalar modes is a necessary consequence of the one-loop effective 
action and trace anomaly of quantum conformal matter or radiation. Thus the conformal anomaly 
provides scalar degrees of freedom in cosmology without the ad hoc introduction of an inflaton. 

These scalar cosmological horizon modes satisfy the full linear response equations to leading order 
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in G N H 2 . In particular, the non-local term is always negligible for these modes, provided the semi- 
classical evolution is begun with modes whose physical wavelength is much larger than the Planck length 



Lpi. Thus the simple analytic solutions to the equations resulting from the anomaly action (7.17) give 
excellent approximations to the exact solutions of the linear response equations around de Sitter space 
to first order in H 2 L 2 pi <C 1. 

Although it is not the full quantum effective action obtained by integrating out the conformal matter 
or radiation fields exactly, the difference is a non-local term which is important for the Planck scale 
solutions of the first kind only. Ignoring this non-local term, the auxiliary field action gives the correct 
linear response equations valid on low energy scales far removed from the UV Planck scale, i.e. in the 
infrared limit where the semi-classical theory can be trusted, and for the second class of state dependent 
solutions. This is consistent with the general arguments that the anomaly action should capture the 
main infrared or macroscopic quantum effects discussed in earlier work |16[ 122] . The fact that there 
is agreement with the local part of the perturbed stress-energy tensor and the infrared relevant terms 
in a context where the results could be compared to exact calculations from entirely different methods 
provides a good indication of the usefulness of the anomaly action as a tool for investigating infrared 
quantum effects due to conformally invariant fields in other contexts where exact results may not be 
readily available. 

In the FRW coordinates the cosmological horizon modes for fixed comoving wavenumber k redshift 
away at late times like a -4 . The linear homogeneous equation these modes obey can also be considered 



in the static coordinates (8.1) of de Sitter spacetime. By recognizing the connection between these 
state dependent solutions to the linear response equations and the solutions to the same auxiliary field 
equations found from the background field analysis in de Sitter space in static coordinates in Ref. [16] , 
we have found that the cosmological horizon modes can describe state dependent changes on the horizon 
scale with a large or even divergent stress tensor at r = H^ 1 (with respect to an arbitrarily chosen 
origin in de Sitter space). This may be understood as a thermal fluctuation of temperature of the 
underlying conformal quantum field away from its value in the Bunch-Davies state of T H = H/2ir. If 
the temperature differs from T H by even a small amount, the corresponding stress tensor diverges on 
the horizon [36J. 

A diverging stress tensor on the cosmological horizon signals the breakdown of the linear approxi- 
mation and the necessity to consider non-linear backreaction effects in the vicinity of the horizon. Thus 
our analysis of the cosmological horizon modes and associated stress tensor perturbations in the static 
coordinates suggests these infrared modes allow the temperature of the causal region interior to the 
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de Sitter horizon to fluctuate away from the BD value, leading to large non-FRW fluctuations of the 
geometry in the vicinity of the cosmological horizon. 
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APPENDIX A: THE LOGARITHMIC DISTRIBUTION IN FLAT SPACE 



There are several ways to work with the distribution K defined by (3.5 ) which lead to similar results. 



The distribution has the same form in a FRW background in terms of the conformal time r] as it does 
in Minkowski spacetime in terms of the usual time coordinate t. In terms of the latter the distribution 
is 



K(t - 1'; k;n) 



dui 
2tt 



-tw(t-f) ln 



-uj 2 + k 2 — lesgnu 
J 2 



(Al) 



As it stands this distribution is well-defined for t ^ t' , but is undefined for t = t'. Two different 
approaches may be considered. In one approach, which owes its origins to Hadamard |38[ [39] . an 



appropriate strictly local distribution is added to (Al) to extend the definition to functions with non- 



vanishing support at t = t' . We shall also describe a second approach, more akin to the Pauli-Villars 



regulator method, which modifies the distribution (Al) at large frequencies and allows for a smooth 
limit at t = t' . 

To begin, a mass m can be introduced into the argument of the logarithm in (Al) by replacing k 2 
by oof. = k 2 + m 2 . Then differentiating with respect to m 2 gives: 







dm 2 



K(t - t';u) k ;n) 



2^ 



-iw(t-f )_ 



1 



sin [oj k (t -t')] 



2 + k 2 + m 2 — ie sgn 



-u 



9{t - t') . (A2) 
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This may be recognized as the spatial Fourier transform of the usual retarded Green's function for a 
scalar field with mass m. For tf > t the retarded Green's function vanishes and therefore is explicitly 
causal. For t = t! the distribution is also defined and vanishes at that point. Note also that after 



differentiation with respect to m 2 , there is no longer any \i dependence in ( A2 ) . The reason for this is 



that the /i dependence of (Al) enters only through the purely local contribution, 

r0 ° du 



V-j- K(t - t';u k ;(j,) 
d\i 



-iuj(t-t') 



TT 



-2 6(t-t'), 



which is independent of k or m. 



If (A2) is integrated with respect to dm 2 = dui 2 , = 2uj k duj k at fixed k, then one obtains 



K(t - t'; uj k ] m) = K bare {t - t'; u k ) + K loca i(t - t'; /i) 



(A3) 



(A4) 



where 



Kbare{t — t;Uf.) 



t-t' 



COS 



(A5) 



This determines the causal K distribution for all t' < t, up to an arbitrary local contribution Ki oca i(t 



t'; fi), which is independent of uj k and satisfies (A3). In the form (A4) one can set m = and oj k = k to 



recover the distribution (3.5). 



The need for the local contribution Ki oca i in ( A4 ) appears when K is integrated against smooth test 



functions f(t') which are non-vanishing at tf = t, for in that case, 



^ dt'K bare (t - t'; k)f(t>) = -2 dt' COS[ ^ / 0] /(^) 

-CXD J —CO 1 ! 



(A6) 



diverges logarithmically at its upper endpoint as £ — > 0. This logarithmic divergence can be removed if 



the [i dependent local distribution satisfying (A3) is taken to be 



Ki ocal (t - t'; n) = -2[ln( M £) + C] S(t - t') 



(A7) 



and the limit £ — > is taken after summing K bare + if/ oca /. Here C = 0.577215... is Euler's constant 
which is a finite additional term added for reason of normalization that will become clear shortly. Then 
one can define the action of the distribution K by the one parameter sequence, 



K{t - t'; k- fj) ~ -2 j [ln(/i£) + C) 5(t - t') + cos [k(t - t')] 9{t - t' - f ) ) 



(A8) 



as £ — > 0. The meaning of the symbol ~ in (A8 ) is that acting upon an arbitrary smooth function f(t'), 



t ( f t_€ dt' 

dt' K(t - t'- k; n) f{t') = -2 lim { [ln(^) + C] f(t) + / — cos [k(t - t')] f{t') . (A9) 



f->o 



t-t' 
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that is, the limit £ — > is to be taken after the integral over t! of a member of a class of suitably smooth, 



test functions f{t') has been computed, provided the integral in (A9) converges 



We note that the definition (A9) is not unique since an arbitrary local distribution Ki oca \{t — t'; fx) 



with the correct dimensions and the property (A3) could have been added to Kb are . The definition flA9 



is in fact an adaptation of Hadamard's Partie finie definition of distributions of this kind [38]. Since 



d dt' 
d£ t - t 



7 cos [fe(t-0] /(0 



cos(/c£) 



fit - 



/(*) 



+ f(t) + o(0 



(A10) 



for test functions which are continuous and differentiable at t = t, it is clear that the integral's loga- 



rithmic dependence on £ is just cancelled by the hi(/i£) term, and the limit £ — > of the sum in (A9) is 
finite. 



This finite limit in (A9) can be demonstrated directly by making use of the cosine integral function, 

r 2n 



ci(z) 



dx 



cos x 



X 



C + ln * + £(")" 



n=l 



2n • (2n)! 



z > 0. 



(All) 



to integrate (A9) by parts. Since 



we have 



d . r . -, cos [w fe (t - 0] 

— Cl \UJk(t — t ) \ = ; 

dt' L v n t-t' 



dt' K(t-t';k;n) f(t') 

o 

-21im\[H^) + C]f(t)-ti(kOf(t-0+ f d£ci[k{t-^)] ^ 



(A12) 



= -2 dt' ci [jfe(i - If)] fj + 2 /(t) In 
provided that /(t') is differentiable, and 



(A13) 



; lim {ci [k{t-t')] f{t')} =0. 



The definition (A9) of the finite part of the logarithmically divergent distribution is closely related to 



dimensional regularization and the subtraction of simple poles at n = 4, corresponding to logarithmic 
divergences of Feynman integrals. As in dimensional regularization and subtraction, the introduction 
of some large but otherwise arbitrary mass scale \i is required for dimensional reasons. The remaining 
dependence on the arbitrary renormalization scale /x in the full linear response equations is removed 



by (2.8), which amounts to the freedom to shift the /i dependence into the unknown coefficient of the 
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local Weyl squared term in the full effective action and the coefficient of the local tensor A a \, in the 
semiclassical linear response equations. 



It is instructive to apply (A9) or (A13) to a simple example of a suitable smooth test function, 
namely, 



f(t) = e^, 7 >0, 



for which the integral (A9) converges, and may be computed analytically. Indeed, 



rt ~t dt' 

2 i — cos[ Wfc (t - t')\ e* = e^Ei [-(7 + + ^ Ei [-(7 - i»k)S] , 



(A14) 



(A15) 



where the exponential integral function for complex arguments, Ei(z) is defined by the analytic contin- 
uation of 



Ei(z) = f 

J — ! 



dx 



x 



n=l 



n ■ n 



(A16) 



for negative real z. Using this in (A15) and (A9) for £ — » gives 

*t 



j dt' K{t - t'; k) /i) e 7 *' = e 7 * In ^ 7 \ 



(A17) 



in which the dependence on Euler's constant has cancelled. The argument of the logarithm is exactly 
K 2 //jl 2 for the function (A14), in which the frequency co 2 is replaced by — 7 2 . This is the reason for 



the introduction of Euler's constant into the definition (A8). In terms of a Fourier transform in time, 



this just amounts to the requirement that the \i introduced in (A7) be identical to the [i in the Fourier 



transform (3.5). The fact that (A17) is strictly proportional to f(t) is a result of the simple Fourier 



transform of (A14), and will not hold in the case of general f(t). 



Note that if 7 > 0, (A17) possesses a smooth limit as k — » 0, and indeed the finiteness of this limit 



can be demonstrated from (A13) for a class of test functions, satisfying 



; hm {ln[n(t-t')]f(t')}=0. 



t'-*-oo 



(A18) 



which excludes fit) = const., i.e. 7 = 0, in which case (A17) diverges logarithmically for k = 0. By 



repeating the steps leading from (A9) to (A13) one obtains for k = 



f dt' Kit -t';k = 0; //) f(t') = -2 f dt' In [/x(t - t')} ^ - 2Cf(t) 
= -2 f dt'\n[ue c {t-t')] 4L, (A19) 



for differentiable functions satisfying (A18). 
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In this treatment of the logarithmic distribution K it has been assumed that the lower limit of the 
integration over t' has been extended to — oo. However, for the application to linear response with a 
well defined initial perturbation at a given finite time to, this definition requires modification. The need 



for a modification can be seen from the integration by parts needed to pass from (A9) to (A13). If the 
lower limit of the integral is replaced by to, there will be an additional contribution proportional to 
ci [k{t — to)] f{to)- From the definition (All) this vanishes as to — > — oo provided that /(io)Ao vanishes 
in this limit. However with to fixed and finite, this surface contribution diverges logarithmically as 
t^t . 

This problem with finite to and indeed the previous divergence as t' — ► t are encountered because 
the non-local kernel has been treated in isolation from the other components of the linear response 
equation, which must be solved self-consistently. In the self-consistent solution the Fourier transform 
has zero support for |u;| — > oo and the full linear response equation contains no divergences for finite 
initial conditions. Since the very large \u\ components will drop out in any case, one may consider the 
associated distribution, 



K PV {t-t';uj k -M) = 



du! 
2?r 



-iw(t-f) ln 



-uj 2 + k 2 + m 2 — iesgnu; 
2 + k 2 + M 2 — ie sgn u 



-uj 



(A20) 



in which the contributions at large |u;| are subtracted against an equivalent distribution for a field with 



large Pauli-Villars mass M. Comparison with (A8) shows that when M S> u) k and K is integrated 
against smooth test functions f{t') with bounded Fourier components possessing vanishing support for 



uj M, the result should become indistinguishable from that of the original distribution in (2.19a), 



with M playing the role of /j, in ( A8 ) . 



The advantage of the explicitly subtracted Pauli-Villars definition (A20) is that its ultra high fre- 



quency components are suppressed, with (A20) vanishing for \u\ 3> M. This leads to the integral 



being well defined for all t — t' with no purely local 6{t — t') contribution of the form (A3). Indeed the 



integration of ( A2 ) between m and M now gives 

2 



K PV (t-t';oj k ;M) 



t - t 



7 {cos[oj k (t - t')} - cosff^O - if)]} 9{t - if) 



(A21) 



with S\ = \Jk 2 + M 2 , instead of (A4). Unlike (A4) this function is strictly non-zero only for if < t, 



vanishing identically at t' = t. It is clear that if M is large, the second cosine in (A21) oscillates very 
rapidly and can give a significant contribution only in the short time interval t — tf < 1/M. Thus, the 
Pauli-Villars regulator plays a significant role only in this very narrow interval of time where if — > t, 
making Kpy essentially equivalent to K\, are outside the region t — t' < 1/M. Although there is no 
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strictly local S(t — t') contribution in the Pauli-Villars method, Kpy has a quasi-local contribution in 



the interval t — t' < 1/M. In this narrow interval the role of the second contribution in (A21 ) is crucial 
since it removes the logarithmic divergence of Kb are , and the difference, Kpy vanishes identically at 
t = t' . Hence the integral over test functions f(t') with non- vanishing support at i' = t becomes 
well-defined. 

For Kpy there is no problem with an initial value formulation and we may define the action of the 
retarded Pauli-Villars kernel on test functions in the finite interval [to , t] by 

«t 



Kpy of = [ dt' Kp V (t-t';k;M)f(t') 

Jtn 



to 



dt' 
t^t 



- {cos[k(t - t')} - cos[ft fc (t - t')}} f{t') 



dt' 
i^t 



7 {cos[k(t - t')] - cos[ft fc (t - t'))} f(t') 



+2 



t0 dt' 



t - t 



7 {cos[k{t - t')\ - cos[O fc (t - t')}} f{t') 



(A22) 



which is completely finite at both its upper and lower limits for suitably smooth test functions f{t). In 
the last form one may recognize the to dependence as giving rise to a temporal transient that falls off 
for k(t - t ) > 1. 



In order to compare the Pauli-Villars definition of the integral kernel (A21) with the finite part 



prescription (A9), (A21 ) may be applied to the same test function (A14). Taking first the infinite range 
(to — > — oo), one finds 

7 2 + k 2 



f dt' K PV (t - t'; k; M)e 7 *' = e 7 * In 



7 2 + k 2 + M 2 



(A23) 



If M 2 3> 7 2 + k 2 this result can be expanded as 

y' dt' K PV (t - t'; k; M)e 7i ' = e' 1 jln (^j^J 



^ + k 2 \ l/^ + k 2 \ 2 



M 2 



+ 



M 2 



+ 



(A24) 



Thus in the limit of large M the Pauli-Villars regularization of the non-local term becomes equivalent to 



the definition (A9), with the identification, M — > fj,, at least as far as the logarithmic dependence upon 



M is concerned. This shows that the arbitrary mass scale of (A9) is equivalent to an ultraviolet cutoff 



scale M, and the additional terms of the expansion in (A24) may be recognized in position space as 
an expansion of higher derivative local terms (— □ /M 2 ) n in ascending inverse powers of the UV cutoff 
M 2 , typical of a low energy expansion of an effective field theory. These finite terms may be neglected 



for sufficiently large M , and are absent entirely in the Hadamard finite part regularization (A9). 
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For finite to the integral formula, 



2 / cos 



[uj k (t-t')} e* = -eT*{Ei[-( 7 + ««; fc )(t-to)] + Ei[ - (7 -iw fc )(t- t )]}, (A25) 



can be used in (A22) to obtain 



jf ^' ifpy (t - t'; fc; M) e 7 *' = |ln 

-Ei[-( 7 + ire)(t-t )] - Ei [_( 7 -ifc)(t-t )] +2ci[M(t-t )]| + •• 



(A26) 



where all terms which vanish for M (7, re) have been dropped. This expression is completely finite 
as t —* to, and in fact vanishes in that limit. Hence there are no unphysical divergences at the arbitrary 



initial time to- This is due to the fact that unlike (Al ) which requires a careful limiting procedure (A9) 



to remove its short time divergence as t' — ► t, the very high frequency response of (A20) is suppressed 



at all times, c.f. (A23) for k — * 00. Moreover, since Ei(z) — > e z /z for large \z\, the last three initial 



state dependent terms in the bracket of (A26) become 



2 e -T(*-*o) [7 cos (k(t - t )) - A: sin (k(t - t ))] sin (M(t-t )) 



+ 2 



(A27) 



t - 1 7 2 + k 2 ' M(t - 1 ) 

The first term decays exponentially for t — to S> 7~ , while the latter falls only as a power times a 
rapidly oscillating function. Thus, at late times, 

ere j\ pyi^i — i ; re; ivi j e ■ — ► e ■ m | 

'to 



dt' Kp V {t - t'; re; M) e 7 *' -> e 7 * In {^^- j , for M » 7, fe , 



(A28) 



which agrees with (A17). Even for constant functions with 7 = 0, the initial state dependent transient 



terms in (A27) fall off linearly with large t — to, and (A28) approaches the Hadamard form (A17) at 



late times, albeit more slowly. This example shows that up to terms that vanish for large M, at late 
times or when operating on continuous and differentiable functions f(t) with bounded first derivative at 
all times, the Pauli-Villars regularization of the logarithmic kernel becomes identical to the Hadamard 
Partie finie definition, with M <-> [x. 

In a different prescription, one could also consider the Hadamard definition of the logarithmic kernel 



for finite to by replacing the lower time limit by to, after the integration by parts in (A13), i.e. 
,-t 



dt' K(t-t';k;n) f{t') 



to 



2 f dt' ci [k{t - t')] fj + 2 fit) In ( k ) . (A29) 

Jto dt 



Computing this for the same test function (A14) gives 
r-t 



£ dt'Kit-t';k; f i)e^ t ' = e 7 * jln 

_ Ei [-(7 + ik)(t - t )] - Ei [-(7 - ik)(t - to)] + 2 e ~ 7 (*-*o) d[k(t - t )] j • 



(A30) 
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This is similar to the corresponding Pauli-Villars form (A26), with \x replacing M and no 1/M terms, 



and is also finite as t — > to- It differs from (A26) only in the last transient term, which falls exponentially 



for 7 > in (A30) but is given by the last term of (A27) in the Pauli-Villars method. 



Unlike (A17) for k = and 7 = 0, the result (A26) remains finite, becoming 

f dt' K PV (t-t';k = 0;M) = -2 In [M(t - to)] + 2 ci [M{t - to)] -2C . 



(A31) 



For times t— to 3> 1/M, the logarithm dominates and (A31 ) grows without bound. Hence the logarithmic 



divergence of (A9) or (A31) for k = 7 = at all times is replaced by a logarithmic growth in time for 



the Pauli-Villars definition of the distribution (A20) 



To see how the logarithmic kernel behaves in position space, one can return to the unregulated form 



of the distribution ( A5 ) which is valid either if t' < t, or if K is integrated against test functions which 



vanish at t = t'. Its form in position space is found to be 



Kbareit) x) 



d 3 k 



Ak-x 



9(t) cos(kt) 



(A32) 



(2tt) 3 [ t 

The integral over k can be performed by going to polar coordinates, and integrating over the angles, 
with the result, 



Kbareiti x) 



29(t) f°° k 2 dk fe ikT -e- il r 



(27r) 2 t Jo 
1 6{t) d 
4ir 2 tr dr j 
1 9{t) d 



ikr 



cos(kt) 



dk 



Jk(t+r) 



+ e 



-ikit+r) _|_ r) 



+ e 



-ik(t—r) 



[8(t + r) + 5(t - r)] , t > 0, r > . 



(A33) 



2-7T tr dr 

Because of 0(t), the argument of the first delta function cannot vanish, except possibly at t = and 
r = 0, where the bare distribution K^ are is undefined. Hence the first delta function in (A33) can be 
dropped if one integrates against test functions which vanish at t = tf and x = x' . 
Reinstating t — t' and x — x' , the effect of the distribution on a test function, 

K bare of = J dt' f d 3 x' K bare (t -t',x- x') f(t', x') . (A34) 

can be computed by first going to polar coordinates centered at x, 

(A35) 



x 



x + nr , 

for some unit vector n, and then introducing the standard retarded and advanced time coordinates with 
respect to (t, x), 

u = t' — t — r , (A36a) 
v = t' -t + r. (A36b) 
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In these coordinates 



Kbare ° / 



oo ;-oo 

d£ / 

-oo J 



r2dr ! ^ojt-t*! a t 



4tt 
dtl 



— e roc 

du I dv 

■oo J —oo 

d 



2vr (t-t')r dr 
d5{v) 



-LlT 
d&n 

4tt 



du 



— oo 



dv 



u — V 

U + V 

u — V 



dv 



f 



U + V 



f 



du 



2 l + d l 
u dv 



v=0 



v=0 



(A37) 



The integral over J_ *\ ^/ diverges logarithmically as e — > + . Thus it may be handled by the same 
Partie finie prescription as in (A9). In order to determine the local contribution at t' = t, x — a 



similar method can be used as that leading to (A19) for fixed k = 0. As in (A18) note that one may 



integrate (A37) by parts and drop the lower limit surface term for functions / for which 

lim {hi(— fj,u)f(u, v = 0)} = , 



u— >— oo 



(A38) 



which excludes spacetime constant functions. Then we may take the limit e — > + and obtain the finite 
result [I0|, 



Kof 



d£ln 

~2tT 



du 



In -n/A )^~ + ~ 
ou 2 dv 



(A39) 



v=0 



for some constant A. Taking into account a relative normalization factor of —l/2ir, (A39) agrees with 



Eq. (20) of Ref. |8j, where the notation H\ was used for the distribution K. See also Ref.|41 
In ( A39| ), / is to be viewed as a function of u, v and h according to 



f{t ,x) = f [t-\ — ,x + n— — 



If one takes the previous example of 



f(t',x') =e^' = e 7*+7(«+«)/2 } 



(A40) 



(A41) 



and substitutes this into (A39), using dxe x hix = —C, one finds 2e 7 ' ln(7//x), in agreement with 



(A17) for k = 0, provided 



1 
A 



(A42) 



Hence, 



Kof = j dt' j d 3 x K(t -t';x- x; n)f(t', x) 



47T 



du 



hi , /xV^VA df + df 

4e / du dv 



(A43) 



v=0 
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for differentiable functions satisfying (A38). 



Unlike (A9), the definition (A39) or (A43) is Lorentz invariant. In the Pauli-Villars regulated 



distribution (A20), which is also Lorentz invariant, the time limits are arbitrary and the lower limit can 



be taken to be finite. Hence condition (A38) may be relaxed. In the definition (A43) there is no a priori 



justification for replacing the lower limit with a finite time boundary. Hence, although equivalent when 
evaluated on functions possessing vanishing support as \uj\ — > oo, which is the physically relevant case, 



the Pauli-Villars subtracted definition (A20) is both well defined over a wider class of test functions, 



and better suited to a finite initial value formulation of linear response. 



APPENDIX B: THE VARIATIONS IN DE SITTER SPACE 



In a conformally flat spacetime with C a b d = 0, and in the conformal "vacuum" state for which the 



background auxiliary fields have the values ip = 0, (p = 2 In $7, the variation of the F a b tensor (4.5) is 
given by: 



SF% = -2(V< a £) (V 6) □ V) - 2(V(» (V 6) U(p) + 2V C [(V C ^)(V°V 6 £) + (V c </3)(V a V fe V)] 

4 
3 



| v a v b [(v c ^)(v c vO] + \ iA(v c ^)(v c v) - 4i? c(a [(v^)(VcV) + (v 6) vo(v c $0] 
^(v a)V 9)(v fe) v) - n(n£>(nV0 + [(v c ^)(v c ^)] + 28\ (V - | 9 crf ) (v c ^)(v^) 



4V c V d (5C {ac b U) -2R cd ^6C a c b d - ^V a V b niP - 4R c( - a (V b) V c iP) + ^R\U^ + | i?V a V 6 ^ 



(Bl) 



Specializing to de Sitter spacetime and using the form of (p = 2Ht for the BD state in de Sitter space, 
one finds 

R 



U(p = -6£r 



2 ' 



(V<^ = 5 a6 (V a ^)(V^) 



-4H Z 



R 



(B2a) 
(B2b) 



Using (B2), the equation of motion (5.14b) for tp and the condition (3.11) 5R = one finds that (Bl) 



becomes 



6F% = -2(V< a £) (V 6) Uil)) + 2V C [(V c V)(V a V^) + (V c ^)(V a V fc V)] - » vav & [(V C ^)(V C ^)] 



» na [(vV)(v c v)] + 1 n (v c ^)(v c vO - ^ (v<°$>) (v 6) ^) + ^5%nip 



-4V c V d (<5C (ac fe) d $ - 2 - v a v b □ ^ + 1 v a v 6 v • 



(B3) 
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Evaluating the a = b = t component of this tensor in the fiat FRW coordinates (3.4) of de Sitter space, 



using (2.2a) and (5.10) gives 



5F\ 



2Ht 5A l 



2 V 2 
3^? 



,2 ■ 



df + Hdt 



2H 2 



AHt (- □ + 4H 2 ) 5R\ =- V 2 v 

3 a 2 



(B4) 



which is (5.23) of the text. 

The variation of the E a b tensor is more involved. Omitting the detailed steps we find in the gauge 



(5.17), 



5E l 



3 * 3 a 2 



d 2 + H& V 



* 2 



2H 2 h tt 



20H Z 



2H 2 -» 
8R\ -—^V 2 {w- 2hu) 
3 a z 



(B5) 



Here w and v are defined in Eq. (5.21). The algebraic manipulation programs Mathematica and 



MathTensor were used to help obtain this result. 
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